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0.1 Abstract 



A phase space theory approach for treating dynamical behaviour of Bose-Einstein 
condensates apphcable to situations such as interferometry with BEC in time- 
dependent double well potentials is presented. Time-dependent mode functions 
arc used, chosen so that one, two,., highly occupied modes describe well the 
physics of interacting condensate bosons in time dependent potentials at well 
below the transition temperature. Time dependent mode annihilation, creation 
operators arc represented by time dependent phase variables, but time indepen- 
dent total field annihilation, creation operators are represented by time inde- 
pendent field functions. Two situations are treated, one (mode theory) is where 
specific mode annihilation, creation operators and their related phase variables 
and distribution functions are dealt with, the other (field theory) is where only 
field creation, annihilation operators and their related field functions and dis- 
tribution functionals arc involved. The field theory treatment is more suitable 
when large boson numbers are involved. The paper focuses on the hybrid ap- 
proach, where the modes are divided up between condensate (highly occupied) 
modes and non-condcnsatc (sparsely occupied) modes. It is found that there are 
extra terms in the Ito stochastic equations both for the stochastic phases and 
stochastic fields, involving coupling coefficients defined via overlap integrals be- 
tween mode functions and their time derivatives. For the hybrid approach both 
the Fokker-Planck and functional Fokker-Planck equations differ from those de- 
rived via the correspondence rules, the drift vectors are unchanged but the 
diffusion matrices contain additional terms involving the coupling coefficients. 

Results are also presented for the combined approach where all the modes are 
treated as one set. Here both the Fokker-Planck and functional Fokker-Planck 
equations are exactly the same as those derived via the correspondence rules. 
However, although the Ito stochastic field equations are also unchanged, the Ito 
equations for the stochastic phases contain an extra classical term involving the 
coupling coefficients. 
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1 Introduction 



Bose-Einstein condensates (BEC) in cold atomic gases are a quantum system 
on a macroscopic scale. At temperatures well below the critical temperature 
almost all the bosonic atoms occupy the same single particle state or mode 
[T]. Consequently, BEC exhibit coherence properties analogous to those for 
an optical laser where macroscopic occupancy of photons in a single mode can 
occur, and therefore their application for interferometry was a logical outcome 
following early experiments demonstrating that BEC could show interference 
effects [7. However, unlike the photons in a laser, the bosonic atoms in a 
BEC interact with each other and this can cause decoherence effects that degrade 
interference patterns. Ideally, the BEC should be at a temperature well below 
the transition temperature, so that coherence is enhanced due to all the bosons 
occupying only one or two modes, with decoherence effects due to thermally 
occupied excited modes becoming small. However such low temperatures are not 
always realised. The interactions between the bosons could in principle be made 
zero via Feshbach resonance methods ^jt^, but this is hard to accomplish over 
the range of trapping magnetic fields involved during an interferometric process 
and requires a broad Feshbach resonance. Furthermore, it is bosonic interactions 
that provide for the possibility of realising Heisenberg limited interferometry via 
spin squeezing effects [5], [6], and proposals for such interferometry have been 
made [7], [8] in which the bosonic interactions are central. Recent reviews of 
BEC interferometry include H], [TO], [U, [13 . 

The theoretical treatment of BEC interferometry must take into account the 
many body nature of the BEC including the presence of boson-boson interac- 
tions. The interferometry process is time dependent and may involve non-static 
trapping potentials, such as when the trap changes from a single well to a possi- 
bly asymmetric double well and back again. The possibility of fragmentation of 
the BEC into more than a single mode must be allowed for, since it cannot 
just be assumed that all bosons stay in one single condensate mode throughout 
the process even if this was initially the case. Decoherence effects associated 
with bosons being lost from condensate modes need to be treated, including the 
creation of Bogoliubov excitations [13] . Unless the BEC is prepared well below 
the transition temperature thermal effects would also need to be taken into ac- 
count. Also, the experiments may not control the total number of bosons in the 
BEC. These considerations suggest that the quantum state of the BEC should 
be described via a density operator rather than a pure state, and the many body 
nature of the BEC would best be discussed in terms of field operators and sec- 
ond quantisation rather than via standard first quantisation methods involving 
symmetrised products of wave functions. 

A number of different theoretical approaches have been used for studying 
BECs, including treatments using variational methods, Heisenberg equations, 
master equations and phase space distributions. General descriptions of and 
references to these methods are set out in [14]. The present paper focuses on 
phase space methods which involve the overall bosonic field rather than just the 
separate modes - these were originally developed for applications in quantum 
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optics [IS] , [H] , [T7] , [IH] , [12] , but are now applied in quantum-atom optics [50] , 
[2T| . The quantum density operator is represented either by a quasi-distribution 
function in a phase space involving variables that replace the mode annihila- 
tion, creation operators for treatments where it is practical to consider separate 
modes, or by a quasi-distribution functional in a phase space involving field 
functions which replace the bosonic field annihilation, creation operators. As 
these two situations are interconvertible it is convenient to consider both though 
ultimately the paper is focused on the latter. Also, it is often convenient to de- 
fine field operators for condensate and non-condensate modes separately [35], 
[23] . and for situations well below the transition temperature often only one or 
two condensate modes are involved. Quantum correlation functions involving 
products of annihilation, creation operators either for the separate modes or for 
the bosonic field are given by phase space ordinary or functional integrals of the 
equivalent products of phase space variables or field functions, weighted by the 
quasi-distribution function or functional. Furthermore, there are several possi- 
ble choices for the type of quasi-distribution function or functional that may be 
used. The positive P type is suitable for consideration of quantum correlation 
functions involving normally ordered products of annihilation, creation opera- 
tors, the Wigner type is better suited to treat those that involve symmetrically 
ordered products. Indeed, when the BEC is well below the transition tempera- 
ture a hybrid approach, such as where the highly occupied condensate modes are 
described via a Wigner quasi-distribution function or functional and the positive 
P type is used for the mainly unoccupied non-condensate modes [20], [24], [25] . 
[14], might best describe the physical situation where the condensate bosons 
behave essentially like a classical mean field and the non-condensate bosons 
exhibit mainly quantum features . In the phase space approaches the Liouville- 
von Neumann or master equation for the density operator is replaced by either 
Fokker-Planck or functional Fokker-Planck equations for the quasi-distribution 
function or functional. Finally, the Fokker-Planck equations are replaced by 
equivalent Ito stochastic equations for either stochastic phase space variables 
or stochastic fields, and the phase space ordinary or functional integrals for the 
quantuim correlation functions given by stochastic averages. 

Phase space methods first involve a consideration of what are suitable single 
particle states or modes that could be used to describe states of the bosonic 
system. In treating both the dynamics of bosonic systems and their static be- 
haviour at non-zero temperatures it is usually convenient to use Fock states as a 
basis. In first quantization the Fock states are symmetrised products where the 
identical bosons occupy single particle states or modes, and for bosons there are 
Fock states in which more than one particle can occupy any particular mode. In 
an idealised Bose-Einstein condensate of non-interacting particles at zero tem- 
perature all N bosons occupy the mode with lowest single particle energy. Finite 
temperature effects for non-interacting bosons allow for excitations to higher en- 
ergy modes. Below a critical temperature Tc Bose-Einstein condensation occurs 
as a phase transition. For non-interacting bosons in a static potential V the 
usual choice for the single particle states is the time independent energy eigen- 
states for a single boson in the static potential. However, for interacting bosons 
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in cases where the trap potential V changes with time - such as in interfer- 
ometry experiments involving Bose-Einstein condensates - there is no obvious 
set of time independent modes that could be used. Of course theoretical treat- 
ments using time independent modes can still be carried out, often these are 
chosen for mathematical convenience such as in facilitating numerical calcula- 
tions. For example Blakie et al [35] use single particle states for non-interacting 
bosons in a static trap potential, Egorov et al [27] use plane wave states. It is 
well known however that for interacting bosons at temperatures close to zero, 
the solution to the time dependent Gross-Pitaevskii equation [25], [H] usually 
provides a good physical choice for the single particle state that all the bosons 
occupy. This mode function is obtained from the Dirac-Frenkel variational prin- 
ciple [3^, [31] in which the state vector is a product with all bosons occupying 
the same single particle state. However, as Leggett has pointed out [1] Bose- 
Einstein condensates may sometimes exist in fragmented states, in which there 
is macroscopic occupancy of more than one mode. Such a situation could occur 
in double- well interferometry with Bose-Einstein condensates, where two modes 
may be involved. These may be two localised modes in the separate potential 
wells, or they may be delocalised modes spread over both wells, in a symmetric 
double well one could be symmetric and the other antisymmetric. Generalised 
coupled Gross-Pitaevskh equations [5 ^ . [55]. [M]. [55]. [55]. [57]. [55 ] . [Tij. [55 ] 
can be obtained via the Dirac-Frenkel variational principle for the actual pair 
of mode functions. Theories in which all the bosons occupy only or a few single 
particle states where the mode functions satisfy Gross-Pitaevskii type equations 
are referred to as mean field theories, since the Gross-Pitaevskii type equations 
contain non-linear terms which can be interpreted as due to mean fields. How- 
ever, the actual dynamics is more complicated than that can be treated via 
mean field theories, processes associated with unoccupied modes such as Bo- 
goliubov excitations [40], [I], [13] need to be considered. Phase space theories 
are a standard approach to moving beyond mean field theories. Since such time 
dependent modes or single particle states provide a useful first approximation 
to treating the dynamics of interacting bosons in time dependent potentials at 
temperatures well below Tc, it is natural to base the phase space theory on time 
dependent modes as these more closely describe the BEG before interactions in- 
volving both occupied and non-occupied modes are included. It is also possible 
that choosing time-dependent modes could facilitate numerical calculations. A 
redevelopment of phase space theory to allow for this modal time dependence 
is warranted. 

However, if the theory is to be based on time dependent modes it follows 
that as the total bosonic field annihilation, creation operators are not time de- 
pendent, then the mode annihilation, creation operators must be. Since in phase 
space theories these mode annihilation, creation operators are represented via 
the phase space variables, the question arises as to whether these should be cho- 
sen to be time dependent to reflect the time dependence of the mode functions. 
It turns out that although the phase space theory can also be developed in terms 
of using time independent phase space variables to represent the time dependent 
mode annihilation, creation operators, the field functions that represent the to- 
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tal bosonic field annihilation, creation operators via mode expansions with the 
phase space variables as coefficients then become time dependent. This is rather 
unsatisfactory as it would be preferable to represent time independent bosonic 
total field operators by time independent field functions. There are other com- 
plications with using time independent phase space variables, so the present 
paper is framed around using phase space variables whose time dependence is 
related to the time dependence of the mode functions. The time dependence of 
the phase space variables is actually chosen to be identical to that for the mode 
annihilation, creation operators. Both the Wigner and positive P distribution 
functions or functionals are treated, both in a double phase space. The theory 
presented here treats two situations, one {mode theory) being where mode an- 
nihilation, creation operators and their related phase variables and distribution 
functions are dealt with specifically, the other (field theory) being where field 
creation, annihilation operators and their related field functions and distribution 
functionals involve a description where individual modes are not distinguished. 

A so-called combined approach could be adopted, in which the modes are not 
divided into condensate and non-condensate types, either when separate modes 
are considered or in the field theory treatment involving the total fields. How- 
ever, even in the combined approach there are some unexpected consequences. 
For the phase space distribution functions involved in the mode theory situ- 
ation, although the Fokker-Planck equations and the results for the quantum 
correlation functions as averages via the phase space integrals are unchanged 
from those for the case of conventional time independent mode functions, the re- 
lationship between the Ito stochastic equations and the Fokker-Planck equations 
now includes an additional classical term that depends on the time dependence 
of the mode functions via the so-called coupling coefficients. These are de- 
fined via integrals of products of mode functions with the mode function time 
derivatives. On the other hand, when the field theory results are derived, not 
only are the functional Fokker-Planck equations for the phase space distribution 
functional and the results for the quantum correlation functions as averages in- 
volving phase space functional integrals unchanged from those for the situation 
of conventional time independent mode functions, but the relationship between 
the Ito stochastic field equations and the functional Fokker-Planck equations no 
longer includes any additional term that depends on the time dependence of the 
mode functions. These results confirm the validity of phase space distribution 
functional results for functional Fokker-Planck equations and Ito stochastic field 
equations presented elsewhere [20] , [2^ in the situation where the total bosonic 
fields are considered. 

The situation changes for the so-called hybrid approach [20], [24], [25], [14], 
which is the approach treated in the present paper. In the hybrid approach, the 
modes are divided into the highly occupied condensate modes and the sparsely 
occupied non-condensate modes, these being treated in terms of a Wigner or 
positive P distribution functions or functionals respectively. The total bosonic 
field operators are the sum of condensate and non-condensate field operators. 
As in Refs [24], [M] time dependent mode functions are involved in both cases, 
and the condensate and non-condensate field operators now also become time 
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dependent. Hence the field functions that represent them in the phase space 
theory are no longer time independent, as was the case for the total fields. The 
phase space theory is therefore more complex than that outlined previously in 
|14j . both when separate modes are considered and when condensate and non- 
condensate fields are involved. However, in both the mode and field theory 
situations, ordinary or functional Fokker-Planck equations and Ito stochastic 
equations for the phase space variables or fields can still be derived, and re- 
lationships between the drift, diffusion terms in the Fokker-Planck equations 
and the classical and noise terms in the Ito stochastic equations established 
both for separate phase variables and for the fields. Hybrid Wigner and posi- 
tive P distribution functions or functionals are treated, both in a double phase 
space. Both the ordinary and functional Fokker-Planck equations now have 
additional diffusion terms involving coupling coefficients between condensate 
and non-condensate modes. The drift terms are unchanged. Also, additional 
classical terms involving the coupling coefficients need to be included in the 
relationship between the Ito stochastic (field) equations and the (functional) 
Fokker-Planck equation, and the noise terms are now related to the new dif- 
fusion term. These additional terms modify the functional Fokker-Planck and 
Ito stochastic field equations presented previously in Ref. [Tl]. We show that 
the condensate and non-condensate stochastic fields involve the same expan- 
sions in terms of time dependent stochastic phases and mode functions as do 
the condensate and non-condensate field functions when expanded in terms of 
non-stochastic phase variables. 

The theory for the hybrid approach is presented in two parts, the first (mode 
theory) dealing with separate time dependent modes, their annihilation and cre- 
ation operators, phase variables, the mode quantum correlation functions, char- 
acteristic and distribution functions, Fokker-Planck and Ito stochastic equations 
for the phase variables associated with separate mode annihilation, creation op- 
erators. This more familiar separate mode theory treatment will be covered in 
Section [5] for the hybrid approach. In the second part (field theory) we con- 
sider field annihilation and creation operators, field functions, field correlation 
functions, characteristic and distribution functionals, functional Fokker-Planck 
and Ito stochastic equations for the fields associated with condensate and non- 
condensate field annihilation, creation operators. This less familiar quantum 
field theory treatment will be covered in Section[3]for the hybrid approach where 
the total fields are divided into condensate and non-condensate fields. However, 
although the quantum field theory treatment is more useful for dealing with 
systems containing a large number of bosons, it is less familiar than the sepa- 
rate mode approach, involving as it does concepts based on functional calculus. 
Consequently the key steps for the quantum field theory treatment will be set 
out as a separate entity, rather than just deriving the functional Fokker-Planck 
and Ito stochastic field equations from the previous separate mode results. Con- 
nections between the separate mode and quantum field theory treatments will 
be frequently emphasised however. Both in the separate mode and the field the- 
ory treatments the stochastic equations are the sum of a classical contribution 
related to the drift vector and a noise contribution involving Gaussian-Markov 
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random noise terms. The stochastic properties of the noise contributions are 
derived and related to the diffusion matrices. The conclusions are presented 
in Section SI along with a comparison to phase space theories based on time 
independent modes. Key results are summarised in Section [SJ both for the hy- 
brid approach and for the combined approach. In Appendix |5] a brief outline of 
functional calculus is presented, and the equivalence of the separate mode and 
field theory treatments is shown in Appendix [7] for the hybrid approach. The 
paper is also accompanied by several further Appendices, available as On-line 
Supplementary Material. Appendices HI M andTUl include details of derivations 
for Sections [5] and [3] to avoid overloading the main body of the paper. Many 
of the technical details on deriving Fokker-Planck equations can also be found 
in Ref. [M]. Time derivatives will be written as ^ though perhaps ^ should 
sometimes be used to emphasise that all time dependences are differentiated. 
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2 Theory - Separate Modes Treatment 



Phase space theory is based on the idea of representing the quantum density 
operator p by a quasi- distribution function P{a, ,a* , which depends on 
phase variables a = {ai, a2, ••, cufc, .., a„}, — {al, aj, .., a^, .., a+}. Here ak 
is associated with the annihilation operator and is associated with the 
creation operator a\ for the fcth mode involved in an expansion of the field op- 
erators ^(r) and ^^(r) in terms of orthonormal mode functions (j)f.{j^,t) or their 
complex conjugates (/<^(r,i). Quantum correlation functions can be expressed 
in terms of phase space averages involving the quasi-distribution function. The 
time dependence of the distribution function is given by a Fokker-Planck equa- 
tion which includes drift and diffusion terms associated with first and second 
order derivatives of the distribution function with respect to the phase vari- 
ables. Ultimately the phase space variables will be replaced by stochastic vari- 
ables satisfying Ito stochastic equations that are equivalent to the Fokker-Planck 
equation, and the quantum correlation determined via stochastic averages. In 
the present paper dealing with time dependent mode functions, a hybrid distri- 
bution function is used where the condensate modes are treated via a Wigner 
distribution function and the non-condensate modes via a positive P distribu- 
tion function. Thus the phase variables for the condensate modes will be treated 
differently to those for the non-condensate modes. The treatment involves using 
time dependent phase variables, and leads to non-standard Fokker-Planck and 
Ito equations. Details of derivations are covered in Appendix [H 



2.1 Field Operators, Hamiltonians, Quantum Correlation 
Functions 

In the Schrodinger picture the field operators associated with annihilating or 
creating a boson at position r are time independent and satisfy the standard 
bosonic commutation rules, for which the non-zero expressions are 

Mr),^\s)]^6{r~s) (1) 

For simplicity we consider only single component bosonic systems, however cases 
where the bosons may involve differing internal (hyperfine) states can be treated 
via appropriate generalisations of the present approach. 

The field operators are fundamental, appearing in expressions for the system 
Hamiltonian 

H = [ dr( — V*(r)1' • V*(r) + *(r)^F*(r) + -^Mt^-M^^M^fr)) (2) 
J 2m 2 

where m is the mass, V is the trapping potential and the zero range approxima- 
tion is used for interactions between the particles, with g — ATrh^ag/m involving 
the s-wave scattering length Os- 

In addition, measurement results are often expressed in terms of quantum 
correlation functions, which involve expectation values of products of field op- 
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erators. A typical quantum correlation function is the normally ordered case 

G"(ri...rp;s,... Si) = (^P(ri)t • • • *(rp)t>i/(s,) • • • *(si)) 
= Tr(*(s,)...*(si)p(t)*(ri)t...#(r,)t) (3) 

where for an N particle system we require p, q < N to give a non-zero resuh. 
Various spatial interference and coherence effects in Bose Einstein condensates 
can be described via such quantum correlation functions. The case for p = q 
where = Sj for all i is proportional to the probability of simultaneously 
detecting bosons at ri , r2 • • • [51] , [3S] . 

2.2 Time Dependent Modes 

For reasons set out in the Introduction, the phase space theory presented here is 
based on using time dependent mode functions to describe single particle states 
of the bosonic system. Such mode functions are required to be orthonormal 
and to form a complete set in the function space of interest. Thus the mode 
functions (j)f.{x,t) satisfy 

J dx(j)l{x,t)(j)i{x,t) = Ski (4) 
^(l)k{x,t)(l)l{y,t) = S{x-y) (5) 

k 

Henceforth x refers to position in a 3D space unless otherwise stated. 

For the case where the set of modes is not complete and restricted to a set 
S the completeness relation is replaced by 

J2M^,t)rk{y,t)^d'ix,y) (6) 

keS 

where S^{x,y) is a restricted delta function [26], [M]. The restricted delta 
function can be used to define a projector that turns any function F{y) into 
a linear combination of the 4'k{x, t) within the restricted function space. Thus 

Vl[F{x)] = j dy5'{x,y)F{y) 

= Y.<PuM I dycj,l{y,t)F{y) (7) 

ktS •' 

If the set of modes is complete then 5^ (x, y) — 5{x — y), the usual Dirac delta 
function. Unless otherwise stated we will consider sets of modes (p^ix^t) or 
4>l.{x,t) that are complete. 
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2.3 Mode Creation and Annihilation Operators 

The field operators may be expanded in terms of a complete set of time depen- 
dent orthonormal single mode functions as 

^{x) =J2Mt)M=o,t), i-Hx) ^Y.^l{t)rk{x.t) (8) 

k k 

where ak{t) ,'al{t) are time dependent mode annihilation, creation operators. 
These operators satisfy standard bosonic commutation rules that are equivalent 
to ([1]) for which the non-zero expressions are 

[ak{t),^l{t)]=Ski (9) 

The annihilation, creation operators can be expressed in terms of the field op- 
erators and mode functions as 

flfcW = dx(t)l{x,t)i'{x) aiit) ^ I dx(t)k{x,t)i>\x) (10) 



so they can be considered as functionals of the mode function or its complex 
conjugate. 

A necessary and sufficient condition for the total field operators to be time 
independent is that the mode annihilation, creation operators evolve as 

^ = Y.Ck^{tra^{t) = Y.ci,mt) (11) 

Ckiit) = J dx^^^^Mx,t) ^ iDkiit) (12) 

showing that the time derivatives of the annihilation, creation operators are just 
linear combinations of these operators. Note that the sum is over all modes. 
By differentiating the orthonormality condition ^ we see that 

Cki + = (13) 

so it follows that the related matrix D is hermitian, D^i = D1^. We will refer 
to the Cki{t) as the coupling coefficients. As will be seen later the coupling 
coefficients play a key role in the functional Fokker-Planck and Ito stochastic 
field equations. 

Assuming that the mode functions 4']^{x, t) or (/'^(a;, t) form a complete set for 
expansions of ^0^(a;, t) or ^0^(a;, t) respectively, we can also expand the mode 
derivatives in terms of the modes with the coupling coefficients as expansion 
coefficients 

d 



dt 

d_ 

dt 



^6l{x,t) = 5]Cfe,</.r(x,t) (14) 
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We can use these results to confirm that the mode expansions for the field 
operators are time independent. We have using (jH)), (|lip and ([T4| 



(15) 
after applying ([T^. 

2.4 Determination of Mode Functions 

There are several possibilities for determining the time dependent modes. One 
possibility is to use the complex conjugates of the eigenfunctions for the first 
order quantum correlation function G^(ri;si) = (^(ri)'l'\I'(si)) - these are the 
so-called natural orbitals. The corresponding eigenvalues give the occupancies 
of these single boson states, and the condensate modes could be chosen as those 
with macroscopic occupancy, with the remainder specifying non-condensate 
modes. However this requires knowing the quantum correlation function, and 
it is hard to see how this could be done without first determining the dynam- 
ical behaviour. Another possibility would be to use the time dependent single 
particle states for a single boson in the trap potential, but this ignores the ef- 
fect of boson-boson interactions. Perhaps the best alternative is to apply the 
Dirac-Frenkel variational principle [30j , |31| to a state in which the bosons are 
restricted to only occupying as few modes as possible, and where the mode 
functions along with amplitudes for the allowed Fock states are treated as vari- 
ational functions. Such an approach leads to the Gross-Pitaevskii equation 
^\ in the case where only one mode is involved or generalised Gross-Pitaevskii 
equations where two or more modes are considered. 

The Dirac-Frenkel principle involves minimising the dynamieal aetion given 

by 

/ m^\^)-{<s>\dt'f)}/2t \ 

S,yn^Jdt\^ -(ci>|iJ|$)/;, j (16) 

where | is the quantum state. For unrestricted states this minimisation 

leads to the time-dependent Schrodinger equation, so the Dirac-Frenkel principle 
enables the fundamental equation for quantum dynamics to be expressed as a 
Principle of Least Action. In the present application however, the quantum 
state is restricted to a specific form involving variational quantities, so applying 
the Dirac-Frenkel principle then does not result in a quantum state that satisfies 
the time-dependent Schrodnger equation, but nevertheless is as close to being 
a solution as that the specific form allows. 

For the two mode case that is relevant to double well BEG interferometry 
the two condensate mode functions are determined as in [M] , [39] . The quantum 
state I of the N boson system as a superposition of the -I- 1 basis states 

I where there are ~ ^ ^^'^ T + ^ bosons (respectively) occupying the 
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two modes with (time dependent) mode functions (f)i{x,t) and (f)2{x,t). The 
amplitude for this basis state is bk(t). 



\m)^ E hit) 

and the basis states are Fock states given by 



(17) 



k)^ \ l' \ |0) ik = ~N/2,~N/2+l,..,+N/2) 



2' / [(f-fc)!]l [(f + fc)!] 

(18) 

These basis states are fragmented or number squeezed states, allowing for both 
modes to have macroscopic occupancy when |fc| <^ N/2. The notation | ^) ^'^^ 
the basis states reflects the feature that the two mode Bose condensate behaves 
like a giant spin system - details are given in [T3] , [35] ■ The total angular momen- 
tum quantum number j = y is macroscopic, and k = — -j,— y+1,--,y~1'~'~T 
specifies the magnetic quantum number as well as 2k determining the difference 
in mode occupancy. 

The quantum state pT|) is a functional of the amplitudes bk and (via (fTO|) ) 
the mode functions (j)i{x,t) and 4>2{x,t). In [14], [39] a self-consistent set of 
equations for the amplitudes and mode functions has been determined by ap- 
plying the Dirac-Frenkel principle [30] , [31] and minimising the dynamic action 
with respect to these variational functions using the state vector given by (fT7|) . 
Here the mode functions and amplitudes act as variational quantities. 

The mode functions satisfy the coupled generalised Gross- Pitaevskii equa- 
tions 
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+ E(5Ei^.mjnC'^J</', (* = 1,2). (19) 

j mn 

and allow for boson-boson interactions and are time-dependent. They follow the 
changes in the time dependent potential V{r, t). This is a generalised mean field 
theory - the quantity in brackets in the final term of the last equation being a 
mean field. The quantities Xij and Yimjn are one-body and two-body correlation 
functions 

X., = (<f|Sla,|<i>) (20) 
Y^^jn = (4>|S|a,|„Sja„ 1$) (21) 

Detailed expressions given in the Appendix for [T3] , [5S] ■, showing that Xij and 
Yimjn are quadratic forms of the amplitudes bk- They are of order N and N'^ 
respectively. 

The amplitudes satisfy coupled matrix mechanics equations 
db 

= E(^fci - ^Uki)bi [k = -N/2, ..,N/2). (22) 
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These A'' + 1 equations ([^ describe the system dynamics as it evolves amongst 
the possible fragmented states. In these equations the matrix elements Hki, 
Uki depend on the mode functions (j)i{X:t)- Detailed expressions for Hki, Uki 
are given in Appendix for [T3], [31]. The quantities Hki are matrix elements 
of the Hamiltonian H in equation ^ between the fragmented states | Y'^)' 
I ^-,1)- The quantities Uki are elements of the so-called rotation matrix, and 
allow for the time dependence of the mode functions. 



Hki 

Uki 



N 



H 



1 

2i 



dt 



N 

y 

N 



N 



N 



N 



(23) 

(24) 



The specific forms of the Xij, Yimjn, Hki, Uki are not important in what follows, 
all that is required is that they have been determined. Equations for the mode 
functions and amplitudes similar to ((T9)) and (|22|) have been obtained by Alon 
et al [37] for single component BECs. The key feature is the self-consistent 
nature of the equations - the amplitudes are required to determine the form 
of the mode equations and the mode functions are required to determine the 
matrix elements in the amplitude equations. 

From the amplitude and mode equations it can be shown that 



ihY,Xij—j dx (f>* {x, t) 4)^ (x, t) 



(25) 
(26) 



The first result shows that the amplitudes would remain normalised to unity 
and the second result is consistent with the modes remaining orthogonal and 
normalised, assuming they were so chosen at t = 0. The second result involves 
the trace of the product of a positive definite invertible matrix X with a matrix 
which is the time derivative of the mode orthogonality matrix. 

For the non-condensate modes a different approach is required since physical 
considerations suggest the mode occupancy would be small. The variational 
approach might be extended to consider quantum states for an iV -I- 1 boson 
system where all except a single boson occupy one or two condensate modes. 
For the case of two condensate modes such a quantum state would be of the 
form 



\m) 



E 



(ai(t)t 



,(f-fc) 



[(f 



k) 



Hi 



[(f 



k) 



111 



{ai{t)^)\Q) (27) 



where the variational functions are now the two condensate mode functions 
4>i{x,t), (j)2{x,t), the non-condensate mode functions (j)i{x,t) and the state am- 
plitudes bk,i (t). Presumably the mode equations for the condensate modes would 
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be similar to those discussed above, whilst the mode equation for the singly oc- 
cupied non-condensate modes would contain mean field terms associated with 
the highly occupied condensate modes - rather like what happens for the high- 
est energy orbital for electrons in an alkali metal atom, the closed shell orbitals 
being analogous to the condensate modes. There would be a multitude of time 
dependent solutions for the singly occupied mode, and these could be used as 
the non-condensate mode functions. Alternatively, the non-condensate modes 
might be just chosen via a Schmidt orthogonality process to be orthogonal to 
the previous condensate modes, starting with a set of plane wave modes or single 
atom trapped modes. 



2.5 Characteristic and Distribution Functions 

Phase space theory involves first introducing characteristic functions that can 
be used to specify all the quantum correlation functions for a given density 
operator. We wish to divide the modes up into two sets - condensate and non- 
condensate modes. The former are to be treated via the Wigner distribution, 
the latter via the positive P distribution. The definition of the characteristic 
function x(^)^^) for this hybrid approach is 

X{^,t) = Tr{n^{^c,^i)^1vci^Nc)P^Nci^Nc)) 

ktNC keNC 

^c{ic,(.c) = expz^(a,^++^,4) (28) 

keC 

where the characteristic variables are ^ = {Ci, •■) ^fe) --^n} = {5c' £,nc}^ = 
{£,1 ^£,2 T-i^k T-^ti} ~ {^C'^jvc }• There are two sub-sets - condensate or 

and non-condensate ^jvc ?jvc 

Using the Baker-Hausdorff theorem 

eyiVi'Y^{akit+£,kal) = (expi^efe4) (expi^afeCfc J (exp(-^ ^Cfe^^) J 

k \ k J \ k J \ k J 

(29) 

and the cyclic properties of the trace, we may combine exponential operators to 
find that 

X{£,e) = exp(-i^e,C+)xp+(€,^+) 

keC 

xp4^,t) = Tr{n+{e)p^-{0) 

n+(C+) = expi J2 n-(0 = expi Yl (30) 

keC,NC keC,NC 

which relates the hybrid characteristic function to that for the case where all 
modes are treated via a normally ordered characteristic function. 
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The characteristic function is related to the quasi- distribution function P{a, a"*", a* , a"*"*) 



via 



exp(i^(afcCfe +^feafe))exp(i ^ a^^^) P(a, a+, a*, a+*) exp(i ^ ^^.0;+) 



fccC 



keNC 



keNC 



d^ad^a'^ exp{i ^ (afe^^ + ^^.a^)) P(a, a*, a^*) 



fc€C,JVC 

where in the second hue all the exponentials have been combined. As in 
[H], [IS], the hybrid quasi-distribution function is of the Wigner type for the 
condensate modes and of the positive P type for the non-condensate modes. 
The phase space variables are a = {ai, a2, ■•, afc, .., a„} = {ac,aNc},ct'^ = 
{af , aj, .., , .., a^} — {aj, a'j^^}, where the set of phase space variables for 
the condensate modes is ac or and that for the non-condensate modes is 



With ckfc = akx + ictky, a]. 



'-kx 



the phase space integra- 



aNC or a JVC- 
tion is 

d^a d^a+ = / / dauxdaky da^^da^y (32) 

k k 

It needs to be emphasised for this double phase space approach that although 
the characteristic function is uniquely determined from the density operator and 
vice versa, the same is not true for the distribution function. As emphasised 
in several references 02], [43] dealing with the positive P representation, the 
distribution is non-unique - many equivalent distribution functions determine 
the same characteristic function and density operator, and these may satisfy dif- 
ferent Fokker-Planck equations based on different choices of the correspondence 
rules used to convert the dynamical equation for the density operator into an 
equivalent Fokker-Planck equation for the distribution function (see Section l^^ . 
Although this non-uniqueness precludes interpreting the distribution function 
as a probability, it is seen as having the advantage of allowing flexibility in the 
Fokker-Planck equations in terms of obtaining related Ito stochastic equations 
that are suitable for numerical work. 

Although the distribution function is not unique, it is known to exist - at 
least for the positive P case. In that case 42 proved the existence of a so-called 
canonical form of the positive P distribution function in the form 



/ 1 \" 1 



a + a 



+* 



a + a 



+* 



B 

(33) 



where 



|a>B = exp(^afcd^) 



(31) 
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are Bargmann coherent states and where the density operator may be written 
as 



a R {a 



+* I 



B \" IB 



Tr{\a)s{a+*\^) 

(34) 



2.6 Quantum Correlation Functions 



The standard approach can be applied to relate the quantum correlation func- 
tions to derivatives of the characteristic functions with respect to i^f. or fol- 
lowed by letting the S,k,^t ~^ ^- This enables the quantum correlation functions 
to be expressed as phase space integrals. The quantum correlation functions are 
for symmetrically ordered products of condensate mode operators and normally 
ordered products for non- condensate mode operators. We have 

G{li,l2,.-lp]mq,..,m2,mi : ji, j2, fcs, ■•, ^2, fci) 

d^ad^a+ P{a,a+,a*,a+*) (35) 

where Zi, Z2, ..Zp; ruq, .., 7712, mi refer to condensate modes and ji, J2, --jr'jks, ■■,k2, fci 
to non-condensate modes. Symmetrically ordered means the average of the 
product of the operators taken in any order. Thus 

" (p + g)! ^ ^^"'1 "'^-"l ""^^-"i)- (^^) 



In Eg. ((36)) the sum over R is over all {p+q)l orderings 5R of the factors --aj o-n 



2.7 Time Dependent Phase Variables 

As outlined in the Introduction, the complex phase variables ak{t), i^) associ- 
ated with the time dependent mode annihilation, creation operators afc(<), aj^(t) 
will be chosen to be time dependent. Hybrid quasi-distribution functions P{a, a+, a*, a+*) 
in double phase spaces will be considered in which the condensate modes will 
be treated in terms of Wigner and non-condensate modes in terms of positive P 
quasi-distribution functions. Note that the distribution function P{a, a^,a* , a+*) 
is not an analytic function of the ak (t) , (t) but also depends on the complex 
conjugates a1{t),a^*(t). The time dependence of the phase variables is arbi- 
trary, but since they are associated with the mode annihilation, creation opera- 
tors one logical choice would be to require the phase variables to have the same 
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time dependence. As we will see, this leads to theoretical results that allow for 
the consequences of the basic time dependence of the mode functions is a rea- 
sonably simple form. In particular, the field functions that are associated with 
the time independent total field operators are also time independent, as might 
be hoped for. An alternative choice of conventional time independent phase 
variables would lead to an unsatisfactory result that time dependent field func- 
tions are associated with the total field operators. The theory can of course also 
be developed via such other choices of time dependence for the phase variables, 
we believe however that the present choice is the most satisfactory. 

The time dependences chosen for the phase variables are thus as for the 
mode annihilation, creation operators 

I I 

and since Cki{t) = iDki{t) with D hermitian, it is possible to find a unitary 
matrix T given by 

^^^=iY.Dkm{t)TUt) (38) 

m 

such that the formal solution for the phase variables is 

au (t) = J2 Tki {t)ai (0) a+ it) = ^ T,* (i)a+ (0) (39) 
( I 

This choice of the time dependence for the phase variables is a key element in 
the derivations. 

By writing ak{t) = akx{t) + iaky{t),a'l{t) = a^^{t) + ialy{t),Tki{t) 
/ [t)+iTj^i (t) in terms of real, imaginary parts, it is then possible to relate phase 
space integrals over the a{t), a^{t) to phase space integrals over the a(0), a"'"(0). 
Thus 

j j d^a{t)d^a+{t) = j j\ldaUt)daky{t)\{da+^{t)da+y{t) 
= J J]ld'^''-(^)daky{0)l[daliO)dal{0) 

k k 

= J J d^a{0)d'a+{0) (40) 

The Jacobian involves the Tki (t) , T^i (t) and the unitarity of T leads to the Ja- 
cobian being unity. Differentiation under the phase space integral can therefore 
be carried out without concern about time dependent differentials - a key sim- 
plification for the derivations. 

The time dependence of the characteristic function variables ^ = {£,i,£,2, --j ?fe 
{it ,^2 ' --'it ^ --^n} is ^1^0 arbitrary. However, if we choose the ^fc(i),C^(i) to 
evolve as for the mode annihilation, creation operators 
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- which is the same as for the phase space variables, then this leads to many 
key expressions involving the characteristic variables being time independent, 
thereby simplifying the theory. 



2.8 Fokker-Planck Equations 

The Fokker-Planck equation gives the time dependence of the distribution func- 
tion P{a, a* , a"*"*). The derivation of it here requires special treatment due 
to the time dependent modes and phase space variables used in the theory. 

2.8.1 Derivation 

We first need to consider the time derivative of the characteristic function (|28p 
in the form given by Eq. ([30]) It follows using Eqs. (fTT|) . (|4T|) together with 
Cki + Cffc = that 



dt 

k 



|E^^«^ = (42) 



9' k 



This means that the time derivative of the normally ordered characteristic 
function Xp+{£,t^~^) is only dependent on the time derivative of the density 
operator, with the time dependences of the mode annihilation, creation op- 
erators ak{t),al.{t) and characteristic function variables C/c(^)j'?fc (0 cancelling 
out. However, although the quantities ^ Ofc^fc a-nd J2 Cfc'^I a-re time 

keC,NC keC.NC 

independent, the same is not true for X^fcec^fc^fc ■ ^^'^ that 

I E ^f^^t -EE (Ci^^^i^t + Cli^k^t) (43) 

fceC keC leNC 

where the /ceC, leC terms cancel out. This means that there are additional terms 
in the equation of motion for the characteristic function. Thus 



keC UNC keC 



keC 

= i\T.T.(^'^^'^l'^t+c*kl^^k^^t))x{^,e)+x{te■,§-/} 

keC UNC 

(44) 

The first term in (j44p then produces extra terms in the Fokker-Planck equation. 
The last term in (|44|) is the characteristic function that would occur if p is 
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replaced by ^pand this term just gives the standard result for the Fokker- 
Planck equation. 

The first term on the right side of (|44l) is given by 



7 UNC 

keC,NC 

X U E E (^'^^ A 7)^ + A 7^ ) ^("' «^ ' ) 1(45) 
X^tciT^c da+ da, da+ dai j 

because using (j3ip and applying integration by parts twice we see that 



d d 

c;2„d2^+ I exp(z V 1 P(a,a+, «*,«+*) 



a, a , a , a 

I rin: nn/.i. 

keC,NC 

(46) 

The second term on the right side of (j44l) is evaluated using the correspon- 
dence rules that apply when -^p evaluated using the Liouville-von Neumann 
or master equation. The derivation of these rules is set out in many textbooks 
(see for example [H], Appendix 12) and here we just present the results. For 
the non- condensate annihilation, creation operators the characteristic function 
is of the positive P normally ordered type so we have 

p ^ akp X(C,C^) =^ TTT^X 
P ^ pa, x(e,e+)^(^ + *e.)x 
P ^ alp x(^,C-)^(^+.a)x 
P ^ pal x(C,e-)=.^X (47) 
For the condensate annihilation, creation operators the characteristic function 
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is of the Wigner symmetrically ordered type so we have 

P ^ pa, + 

P =^ alp X{U^)=^[Q^^+I^^t)x 

P ^ pal X(e,e-)^(^4^e^)x (48) 

From these results the corresponding changes to the distribution functions can 
be deduced using pSI) or ((3T|) . This involves integration by parts. The ana- 
lyticity of the functions expi ^^'^^^{^^.0;^} and exp i ^)!^-^{Q!fc^^} enables two 
options for partial differentiation {atx or iaty and or ia'^y),to be obtained, 
but this need not be made explicit here (see [42], [43] or [14] ). 

The correspondence rules for the non- condensate annihilation, creation op- 
erators are 



P => akP 



P{a, a"*", a*, a^*) ak P 



p ^ pdk P(a,a+,a*,a+*) ^ ( -7;— p + Qffc ) P 



d 



p ^ alp P(a,a+,a*,a+*)^ -— +a+)F 



d_ 

dak 

p =^ pal P{a,a+,a*,a+*)^ alP (49) 

The correspondence rules for the condensate annihilation, creation operators 
are: 

p cikP P(a, a+,a*,a+*) (afc + ) F 

p => pcik P(a,a+,a*,a+*) ( afc - ) P 



2 a. 



p ^ alp P{a,a+,a*,a+*)^[a+--^^^P 

p ^ pal P{a,a+,a*,a+*)^ (^a+ + ^^^P (50) 

where the a, a^, a*, a"*"* are regarded as four independent complex variables. 

By applying the correspondence rules ([49]) or ([50l) in succession to the various 
products of the density operator with mode annihilation, creation operators that 
arise from the term the second term on the right hand side of Eq. (|44)) can 
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be written as 



[ [d'ad'a+eM^ E + W)) j- E(9^^a 

keC,NC V k 

xP{a, a+, a*, a+* 

d^ad^a+ exp{i E K?fc+?fcafe)) 



do^k 



^EE(' 



fceC,AfC 



xP(a,a+,a*,a+*) (51) 
To obtain the specific forms of the drift Af^ , A'j^ and diffusion terms D^^ ^ , D^^ ^ , D^^ + , D^^ ^ 



the Liouville-von Neumann equation 

^h■^^p^[H,p] (52) 

is used in conjunction with the correspondence rules, with the Hamiltonian 
obtained from Eq. ([2]) by expanding the field operators in terms of mode anni- 
hilation, creation operators using Eqs. (jS]). The Hamiltonian will then contain 
terms of the form a\di involving one mode creation operator and one mode an- 
nihilation operator from the kinetic and trap potential energies, as well as terms 
of the form al,ajdmCin involving two mode creation operators and two mode an- 
nihilation operators from the boson-boson interaction energy. These terms will 
be multiplied by spatial integrals involving the mode functions - the latter be- 
ing assumed known (see Section [2T4|) . Division of the modes as condensate and 
non-condensate is needed to identify whether Wigner or positive P correspon- 
dence rules apply. For the condensate modes, third order derivatives occur in 
the Fokker-Planck equation, but these are discarded as being too small to mat- 
ter because of the large occupancy of the condensate modes. It is well known 
that there will only be first and second order derivatives associated with the 
non-condensate modes, these being associated with the positive P distribution. 

We also need to consider the time derivative of the characteristic function in 
the form Eq. ((3T|) . Differentiation under the phase space integral can be carried 
out without any issues associated with time dependent differentials, as Eq. (I40p 
shows. Also we see that similar to the proof of (j42|) 

I E ^^^i -§iT. - 1 E("fc^^ + = (53) 

k k k 

Hence we have 



-1- * \ 
a,a ,a ,a ) 



ktC,NC 



(54) 
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Substituting the results from (gS]) and (HI]) into Eq.(|31]) the Fokker- 

Planck equation for the distribution function can be obtained in the form 



As the expressions are now becoming rather lengthy we now change to a simpler 
notation. 

2.8.2 Notation Change 

We first introduce the symbol A to designate condensate and non-condensate 
modes, where A = C,NC. The mode functions (j)f.{x,t) and their conjugates 
(pKxjt) will now be designated 4)j^j.{x,t) and their conjugates (l)\j^{x,t), where 
fc = 1, 2, .., riA lists the separate condensate or non-condensate modes. Both the 
(j>Ak{^j </'^fc(x, t) may finally be listed as (t>Aki^^ '^^^'^ = — , +, so that 

(j)~^f,{x,t) =(f>j^i,{x,t),(f>\f,{x,t) =(f)\^{x^t). Note that if there are nc condensate 
modes and non-condensate modes, the total number of (j)'^j,{x,t) will be 
2(nc + npfc), which is twice the total number of modes. This is because there 
are n = nc + n^c modes 4>i.(x,t) plus their complex conjugates (j)l{x,t). The 
phase space variables will now be written as a^Afe and those in the characteristic 
function as ^^^fc- The equations for the phase space and the characteristic 
function variables now become 




(55) 




(56) 



Bl 




(57) 



Bl 



with 



C 



AkBl 



^Ak Bl (^) 




(58) 



AkBl 




(59) 
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Note that the sums in and (1571) include both condensate and non-condensate 
modes. The original drift, diffusion terms are now designated as ^Ak^^AkBi- 
The time derivatives of the modes are now written as 

§-M^,t)^j2^Ai:BirBi{^,t) (60) 

Bl 

The orthonormality and completeness relationships will now be 



^(t)Aki^^t)^Akiy^*) = St^Aix,y) (61) 

k 

where 5^a{x, y) is the restricted delta function for the fiA modes and the pro- 
jectors r^^ now are 

K^'lFix)] = I dySMx,y)Fiy) 

k Bl k 

B^A 

SO that the projector V!^^ acting on an arbitrary function written as a linear 
combination of the modes 0^^(a;,t) and orthogonal modes (j)g[{x,t) [B ^ A) 
projects the function onto just the linear combination of the modes (j)'^i^{x,t). 

2.8.3 Fokker-Planck Equation Result 

In terms of the new notation the Fokker-Planck equation will be of the form 

- a^'^Ak + ^ E E a^a^'^Ak bi ] ^(«, a\a*,a^*) 
t^k ^"^^'^ 2 t^k t^i ^"'^■^/c da.Bi J 

(63) 

where A^^, is the drift vector and i?^^. gi is the diffusion matrix. The diffusion 
matrix is given by 

^AkBl = F)^AkBl + ^{^ AC 5 B NC 5 ^ -uCXk Bl + ^ B C 5 A NC -tiC'^^l Ak) (64) 

The first term D'^^^^ is from the standard derivation of the Fokker-Planck 
equation via applying the standard correspondence rules to the terms in the 
characteristic function that arise from ^p, the second arises from the time 



24 



dependence of the condensate and non-condensate mode functions. Note that 
the diffusion matrix is still symmetric 



'Ak Bl 



BlAk 



(65) 



The drift A'^j^ vector and the diffusion matrix E^^^^ are functions of the 
otiiAk- In the present approach using time dependent phase variables the drift 
vector A^j, vector is the same as those that would be obtained if time inde- 
pendent modes were used, but the diffusion matrix E^f, is changed from the 
expected . The Fokker-Planck equation can also be expressed in terms 

of real variables which involve the real and imaginary components of the phase 
variables a^^Afe , but we will not do that here. In deriving the Fokker-Planck equa- 
tion there are often terms involving third and higher order derivatives arising 
from the condensate modes and their Wigner representation. These are usually 
small, scaling as higher powers of 1/Vn, and hence are discarded. Thus we see 
that the Fokker-Planck equation has changed in this hybrid approach from that 
for time independent modes. 

2.9 Langevin (Ito) Stochastic Equations 

In this section we determine the Langevin equations that are equivalent to the 
Fokker-Planck equation for the distribution function. These will be in the form 
of Ito stochastic differential equations. The derivation is based on that given by 
Gardiner [43j . but modified to allow for the phase variables being time depen- 
dent. The final stochastic averages would determine normally ordered quantum 
correlation functions associated with non-condensate operators, or symmetri- 
cally ordered quantum correlation functions associated with condensate oper- 
ators. Note also that the Langevin equation derivation does not depend on 
the distribution function having any particular properties, such as being real or 
positive. We must of course assume it to be non-analytic in general. 

2.9.1 Phase Space Average 

The phase space average of the functions, F(a, a^) is given by 



For determining the quantum correlation function G(Zi,^2, ■■lp]fnq, ■■,m2,mi :: 
Jii j2, ■■jr]ks,..,k2,ki) in Eq. ((35| the function is 

F{a,a+) = [a/^ a™^..am2a„J[a+ afe^..afc2afcj. The phase 

space average will change with time not only because the distribution function 
P{a, , a* , a~^* ,t) is time dependent, but also because F{a,a'^) is now time 
dependent since from Eg. ([57)1 the a,a~^ depend on time. Note the differentials 
J J dPa^iPa are time independent. 




(66) 
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The change in the phase space average from times ttot + 5t due to ah factors 
is given by 

(F(a,a+))_^-(F(a,a+))^ 



d 



d^aF{a,a+) \ —P{a,a+,a*,a+*,t) \ ] St 



d 



d-'a+d-'a ^—F{a,a+)j P{a, a+ , a* , a+* ,t) j St (67) 

The second term in the last equation is absent in the standard treatment [33] 
based on time independent phase variables. This result gives the change in the 
phase space average correct to 0{St). 
Now from ([5^ 



d_ 

dt 



2^ F{a,a+)- 



fiAk 



d 



dt 



J2 J2 ^Ak Bl O^f^Bl 



(68) 



Bl 



Substituting for -^P from the Fokker-Planck equation and for -g^F 
from we get using integration by parts and assuming that the distribution 
function goes to zero fast enough on the phase space boundary 



(F(a,a+))_^-(F(a,a+))^ 



t+st 
d^a+d^a 



E 

fj.Ak 

xP{a, a'^ ,a* ,a'^* ,t)St 



da^j_Ak 



F{a,a+] 



Alk+T.^'lkBl^ 



fj.Bl 



Bl 



Sa^ Sol 



fj,Ak uBl 



d d 
da^Ak da^Bi 



F{a,a+) 



^AkBl 



xP{a,a'^ ,a*,a'^*,t)St 



(69) 



where the terms involving first order derivatives have been combined. 
Hence 



fj.Ak 



-F{a,a+) 



da^Ak 



fj,Ak uBl 



dUf^Ak davBi 



^Ak + E '-^Ak Bl O'f^Bl 
Bl 

F{a,a+) 



^Ak Bl 



(70) 
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giving the time derivative of the phase space average of F(a,a^) as the phase 
space average of the sum of the two quantities in the {} brackets. Note that no 
specific properties of the distribution function were needed. This result will be 
used later to derive results for time derivatives of quantum correlation functions. 

2.9.2 Stochastic Variables Average 

For the Ito stochastic approach we now replace the variables a^Ak by stochastic 
variables a^^^.. In terms of these stochastic variables the stochastic average at 
time t of F{a, a+) is given by 

F{as{t),a^+{t)) = j^Y.fKAkm) (71) 

i=l 

where a^^fci(i) is the ith member of the stochastic ensemble of N samples and 
/(a^^j. (t)) is the same as F{a^{t),a^^{t)). The key idea is that the phase 
space average at any time t of arbitrary functions F{a^a^) and the stochastic 
average of such functions are made to coincide when the stochastic equation for 
the Q:^^fe(i) is suitably related to the Fokker-Planck equation for the distribution 
function P{a,a'^ ,a* ,a'^* ,t). Thus 

(F(a,a+))^ = j j cpQ+(faF{a,a+)P{a,a+,a*,a+*,t) 
= F{a'{t),a'+{t)) 

= JfT^fKAkiit)) (72) 

In turn, the phase space averages are related to quantum averages. In particular 
the normally ordered quantum correlation functions are given by the stochastic 
average of the product of stochastic c-number phase space variables 

«l -"l"-, -am.am,)^ = «,(i)-a^,(i)) «'(*)••<(*)) (73) 

For bosonic systems such stochastic averages involving c-numbers can be carried 
out numerically, and this method is often more efficient than having to determine 
the full distribution function. . 
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Now the difference in the stochastic average from time t to t + St is given by 



F{a^{t + St), a^+{t + St)) - F{a^{t), a^+{t)) 
{F{a^{t + St),a'+{t + St)) - F(a«(i), a«+(i))} 




F{a,a+) 



fxAk vBl 



(74) 



where we have used a Taylor expansion for /(aj^^fc(i + St)) with the notation 



5aUk(t) = al^Akit + St) - a^A^it) 



(75) 



for the fluctuation in a^^^.. In obtaining the last result the stochastic average 
of a sum being equal to the sum of the stochastic averages has been used. It is 
ncccHsary to consider fluctuations of the second order because - as we will see - 
stochastic averages of such terms are of 0{St). 

Now suppose oi^^Aki^) satisfies an Ito stochastic equation of the form 



<^lAk{i + Si) - o^^Akii) 



= KkHDmmt + T.^AkaKDmit)) [ 



t + St 



dtiTaiti) (76) 



which gives the change in Q!*^fe(t) correct to 0{St). The Ito stochastic equation 
is also written as a differential equation 



d 



A'AkKn^it))+22^AkaKD^it))'^a{t+) (77) 



where the stochastic behaviour is due to the Gaussian-Markoff random, noise 
terms Ta- The aim is to find expressions for the A'^^. and the B'^f,^ so that the 
phase space and stochastic averages coincide for an arbitrary choice of F{a, a+). 

The Tait) are real Gaussian-Markoff random noise terms (a = 1, 2, ..) whose 
stochastic averages are given by 



Ta{ti 



ra{ti)rbit2) = 
TJMW^Jt2)T\At3j = 





SabSitl 





t2) 



Taitl)Tb{t2)T^it3)rd{t4) = ra{ti)Ti,{t2)T4t3)Td{ti)+Ta{ti)rc{t3)Th{t2)Tdit4) 

(78) 



+Ta{ti)Ta{t4)n{t2)T^{t3) 
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so that the stochastic averages of products of odd numbers of F are zero and the 
stochastic averages of products of even numbers of F are the sums of products 
of stochastic averages of pairs of F. For the moment we leave the number of Fq 
unspecified, the number will turn out to be 2{nc + unc)- It is also assumed 
that any function h of the a^^j, (t) at time t and the F^ (t) at later times are 
uncorrelated - this is the meaning of the Fa(t+) in (177)1 . Thus 



= Haf^^^ih)) rait2)rb{t3)T,{U)..Tk{tl) h <t2,t3,..,tl (79) 

With these results we can now obtain expressions for the stochastic averages 
in Ea. (l74|) . Details are given in Appendix |8l 

For the first order derivative terms we find that 




(80) 



where the stochastic average rules for sums and products have been used, the 
non-correlation between the averages of functions of a^^^, (t) at time t and the 

F at later times between t to t + 5t is applied, and the term involving Ta{ti) is 
equal to zero from ([78|) . Note that the first order term is proportional to St. 

For the second order derivative terms, we have on expanding the product 
^'^liAki^)^'^tBii^) > using the stochastic average of a sum being the same as the 
sum of stochastic averages, the result that the stochastic averages for the func- 
tions of the ct^u„iit) and the Taitj^.) are uncorrelated, and after using the result 
that terms Fa(ii) involving single F's give zero together with the result for terms 
Fa(ii)Ff,(t2) involving two F's 



/ dt^Taih) / dt2Tt{t2) = dti dt2 Ta{tl)Tb{t2) 

Jt Jt Jt Jt 

/t+5t pt+5t 
dti J dt2 Sab6{ti - t2) 



rt+St 



ft+St 



Sab St 



(81) 
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we find that the second order derivative term is 



nAk uBl 



d d 
da^Ak da„Bi 



Fia,a+) 



p,Ak uBl 




T.^'AkaHD^ii))^maHD^^)) 



St 



^EE 

p,Ak uBl 



d 



d 



da^Ak da^Bi 



F{a,a+) 



St 



(82) 



after neglecting terms of order St'^ Note that the second order term is also 
proportional to St 

The remaining terms give stochastic averages correct to order St^ or higher 
so that we have correct to first order in St 



F{a'{t + St), a^+it + St)) - F(a'*(i), a'+{t)) 



E 

fiAk 



d 



da^j_Ak 



F{a,a+) 



+ UEE 



nAk uBl 



d d 



dafj,Ak da 



uBl 



-F{a,a+) 



malnJt)) S^(a|^,„(t))]^;'r,s, 



St 



(83) 



dt 



F{a'{t),a'+{t)) 



E 

p,Ak 



d 



da^Ak 



F{a,a+) 



■^Aki^lDrnit)) 



1 



^EE 

fiAk uBl 



d d 
_dap,Ak dauBi 



F{a,a+) 



(84) 



This result is exactly the same as for the standard treatment based on time 
independent mode functions, since the Ito stochastic equations (177)1 have exactly 
the same form in that situation. 

The result ([M)) for the stochastic average will be the same as that in (|70l) 
based on the phase space average if we have the following relationships between 
the matrices A and E in the Fokker-Planck equation and the matrices A and B 
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occurring in the Ito stochastic differential equation. 

Bl 

EZbi (85) 

If the modes were time independent, the couphng constants would all be zero 
and the relationship would have been A'^^. = Aj^^, and [B B^]'^^ gi = D^^^i, 
which relates the Ito equation quantities A and B to the drift A and diffusion D 
terms that occur in the standard Fokker-Planck equation obtained from just the 
■^p term in the time derivative of the characteristic function. This is the usual 
condition found in the textbooks |43j . For the hybrid case with time dependent 
modes not only is the relationship changed for the diffusion terms, but it is also 
changed for the drift terms. The difference is entirely due to the presence of 
the coupling terms C^j. that depend on the time dependences of the modes 
(/)^^(a;,t). Clearly, if the drift vector A from the Fokker-Planck equation and 
the matrix C from the time dependence of the mode functions are known then 
the vector A in the Ito stochastic equation can be found. It is known [55] that 
the complex symmetric matrix E may be factorised in the form KK^ = E, so 
the matrix B in the Ito stochastic equation can also be determined. 

B = K 

KK^ = E (86) 

This result is known as the Takagi factorisation 35] • The proof is given in [46] 
(see section 4.3). The construction involves the eigenvectors of the matrix EE*, 
which is hermitian because E is symmetric, and which also has non-negative 
real eigenvalues This result is not well-known, and does not require E to be 
positive semi-definite, as is sometimes thought to be the case. Note that in 
general S is a complex 2n x 2n matrix. Also note that B is not unique, since 
with any orthogonal matrix TZ we also have {BTZ){BTZ)'^ — E. Again, the total 
number of Gaussian-Markoff noise terms Ta(t^) is 2n, the total number of aAk 
and ce~^j., or equivalently twice the total number of modes. 



•^Ak 



[B S^i^'"^ 



iAk.Bl 



2.10 Expressions for Ito Stochastic Equations 



We can also write the Ito stochastic equation (|77|) in terms of quantities that 
appear in the Fokker-Planck equation. From Eq. (|85|) we have for the Ito 
stochastic phase variable equation in the hybrid case 

- ^Ak + E ^Ak Bl <Bl + E ^Aka 
I a 

where A is the drift vector and K is related to the diffusion matrix E via 



J2^AkaK^Bla = 



AkBl 



^Ak BI+ 2 ('^^ Sb NC ^ti -uCXk bi + SbC Sa NC S„ -t^Csi ^fe) 
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This form of the Ito equation is the most useful as it involves the drift, diffusion 
terms in the Fokker-Planck equation ((63)) plus the coupling coefhcients that 
allow for the modes being time dependent. Note that the diffusion terms (|64p 
in the Fokker-Planck equation also involve the coupling constants. Thus the 
equivalence between the quantities A, B in the Ito stochastic equations and 
the drift A, diffusion D terms in the standard Fokker-Planck equation and 
the mode time dependence matrix C can be established. There is an extra 
term C^^ ct^s/ (which is linear in the phase variables) in the drift term, 

and an extra term \ {5ac5bnc (5^ -^^Ak bi+^bc Sa nc -t^C^i Ak) (which is 
independent of the phase variables) in the diffusion term, which are not present 
in the standard treatment involving time independent modes. 



2.11 Classical and Noise Terms 

We can write the Ito stochastic equation in terms of a classical term and a noise 
term 

d f d , \ f d ^ 



Ql-,Ak - \jt^^^^)^^^J\jt^^^'^ 

Ol^tAk j = ^Ak + X] ^Ak Bl ^IlBI 

/ class I 

= E^A.ara(t+) (89) 

/ noise 

If only the classical terms were included, then the solution for a^^;. would 
determine classical trajectories in phase space; hence their name. 



2.12 Properties of Noise 

The stochastic averages of the noise terms can now be evaluated. For a single 
noise term and the product of two noise terms we have using Eqs. (|75|) and 

(^°--"'),. = ° 

(^"M>(i.)) (^"""'"'O = - '2)E5U,('.,2) (90) 

showing that the stochastic average of a single phase noise term is zero whilst 
that for the product of two phase noise terms is delta correlated in time and 
equal to the appropriate diffusion matrix element. 

In fact the diffusion matrix elements determine all the stochastic averages 
of products of phase noise terms. With an odd number of terms the stochastic 
average is zero. For an even number of terms the stochastic average involves 
sums of stochastic averages of products of diffusion matrix elements, which 
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reflects the Gaussian-Markoff properties of the Ta- Thus for three and four 
phase noise terms 



(^"^(ii)) (^<B/(i2)) (|a|c.„fe)) =0 (91) 



{^t^lcM)^ 

+Sih - U)Sih - ia) E^^JhA)E^iZJI^) (92) 
Note that the phase noise terms are not themselves Gaussian-Markoff processes. 
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3 Theory - Quantum Field Treatment 



In this section the field operator and functional phase space approach is de- 
veloped for the situation treated in Section [2l where the modes are divided 
into condensate and non-condensate modes. In general, the phase space func- 
tional approach is based on representing quantum density operator by a quasi- 
distribution functional in a phase space involving field functions which replace 
the bosonic field annihilation, creation operators. Functionals and key proper- 
ties are summarised in Appendix [6l In this hybrid approach time dependent 
condensate and non-condensate field operators are involved and these are repre- 
sented by time dependent condensate and non-condensate field functions. The 
distribution functional will be of the Wigner type for the condensate field and of 
the positive P type for the non-condensate field. As shown in Appendix [71 this 
treatment is equivalent to the treatment in Section[2]based on quasi- distribution 
functions of phase space variables that replace the annihilation, creation opera- 
tors for the separate modes. Details of derivations are contained in Appendices 
[9] and [101 



3.1 Condensate and Non-Condensate Field Operators 

Field operators ^c(a^) and ^nc{x) for condensate and non-condensate modes 
respectively may also be defined in which the sums in ([8]) are over restricted sets 
of modes '14 . The condensate and non- condensate field operators are 



^c{x,t) = Y.Zkmk{x,t), =^St(t)0*.(x,i) 

feeC fceC 

i>Nc{x,t) = ^ = I] (93) 

keNC keNC 

where the mode sums are restricted to be only over condensate or non-condensate 
modes. Their sum gives the total field operators 

i'{x) ^^cix,t) + i'Nc{x,t) i'Hx) ^¥^,{x,t) + ¥j^c{x,t) (94) 

which satisfy ([8|). 

These operators satisfy commutation rules involving restricted delta func- 
tions (see dH])). 

[■i/c{x,t),^l{y,t)] = Scix,y) [^Nc{x,t),i'%ciy,t)] = ^Nc{x,y) 

[i'c{x,t),¥j^ciy't)] = [i'Nc{x,t),i>l,{y,t)]=0 (95) 
Sc{x,y) + Sc{x,y) = d{x-y) (96) 

and condensate and non-condensate operators commute. In terms of projectors 

^c{x) [*(x)] ^Nc{x) ^V^^[i'{x)] (97) 

We note that the condensate and non-condensate field operators are time 
dependent - it is only their sum that is time independent, as shown in ()15p . The 
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feature that the separate field operators are time dependent requires significant 
modifications to the standard derivations of the FFPE and Ito stochastic field 
equations. We find that 

§-^^cix,t) = ^ ^ C,,a,(i)0fc(x,t)+^ ^ aKi)Crfe0fe(x,i) 

keC UNC UC keNC 

+ (98) 

because now there are pairs of fc, I such as fceC, leNC for which the cancellation of 
Cki via the required C*^, using cannot occur, since for the latter leC, keNC. 
The other three condensate and non-condensate field operators are also time 
dependent. 



3.2 Condensate, Non-Condensate and Total Fields 

The total field operators ^(x), ^'^{x) in Eqs. ^ are represented by total field 
functions, which are defined by equations analogous to ([5]) for the field operators 

'4'{x,t) =^ak{t)(pk{x,t) ^i^ix) ^^{x,t) ^^al{t)(j)l{x,t) =-0+(x) 
fe fe 

(99) 

where au (t) , (t) are time dependent mode phase variables that satisfy Eq. 
([57)) . Using Ea.([ni) and Cki + Cl^. = we see that the total field functions are 
actually time independent. 

-i,{x,t)=Q -4,+ [x,t)^0 (100) 

This result depended on all of the (p^ix, t) or <j)*k{x, t) being involved in the field 
functions, since expanding ■^(f>f.{x,t) or ^(/)^(a;,t) involves all modes. Modifica- 
tions are needed when restricted sets of modes are involved for condensate and 
non-condensate field functions. 

However, in the hybrid approach where condensate and non-condensate modes 
are treated differently, the phase space functional approach is based on repre- 
senting quantum density operator by a quasi- distribution functional in a phase 
space involving four field functions which replace the bosonic field annihilation, 
creation operators ^c{x, t), "^cix, t),"^/Nc{x, t), '&]yp(x, t). The condensate and 
non- condensate field functions are defined by equations analogous to Eq. (j93p 
for the field operators 



i^c{x,t) = ^ak{t)(l)k{x,t) ip'^{x,t) ='^a+{t)(f>l{x,t) 

ktC keC 

i^Nc{x,t) = ^ ak{t)(j)k{x,t) ^l^'l^fjix^t) ^ ^ a'l{t)(j)l{x,t){l{)l) 

keNC keNC 

where au {t) , {t) are time dependent mode phase variables that satisfy Eq. 
P7p. In this hybrid case the field functions are time dependent. Similar to Eq. 
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the time dependences of the field functions are 
d 

keCleNC UC kcNC 



dt 



i>c{^,t) = Ckiai{t)(j)k{x,t) + Mt)CtkMx,t) 



keC leNC IcC keNC 



d 

TrM^Nci^^i) = E E'""°''W'^fe(^'*) + E E°''W'^*fc'^fc(^'^) 



keNC UC UNC keC 



dt 



dt 



keNC leC leNC keC 

(102) 



However, in this case there is incomplete cancellation because Cki and C;*^ in- 
volve different pairs of fc, I so these derivatives are non-zero. 

3.3 Characteristic and Distribution Functionals 

The field theory phase space approach involves introducing characteristic func- 
tionals that can be used to specify all the quantum field correlation functions 
for a given density operator. We first define four time dependent characteristic 
field Junctions via 



keC keC 

^Ncix,t) - J2 ikit)M^,t) = E CW<^fc(2^'*)(103) 

keNC keNC 

with the same time dependent c-numbers Cfe(i):Cfc (^) before. These satisfy 
Eq.(j4T|). These sum to total characteristic field functions 



S+=S++S+c E-=E^,+E^C (104) 
The hybrid characteristic functional x[Sc, S^, S^vci S^^] is defined via 

S^c,S+d - Tr{n'^[Ec,E+]n+[E+c]pn-[ENc]) 
^^[^wc] = «^^P* J d.x-^Nc{x,t)Ejj^{x,t) 

Ci~[Enc] = expi J dxENcix:t)ii^j^^{x,t) 

n^[Ec,E+] = expi J dx{^c{x,t)E+{x,t)+Ec{x,t)^l,ix,t)) 

(105) 
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By applying the Baker-Hausdorff theorem we find that 
r2+[Sj] = expi J dxi'cix,t)E^{x,t) 

h'[Ec] = expi J dxEcix,t)i'l.{x,t) (106) 
we find that the overall characteristic functional is given by 

x[Sc,S+,Sa,c,S+ J = exp|-i J dxEc{x,t)E+{x,t)^ Xp+[Sc, S+, Sjvc, S+ J 

(107) 



where 



Xp+[^c, sJ, Enc, =^nc] 



= rr(|expi J dx{^c{x,t)E+{x,t)+^Nc{x,t)E+f.{x,t))^ p 
X jexpj J dx{Ec{x,t)^lj{x,t)+ENc{x,t) ^l^(j{x,t)) 

. rK{»p./<i.,*,.)H.W)}^{exp./..,HW*.M,} ,108) 

is an auxiliary characteristic functional. Similar to the separate mode situa- 
tion, this characteristic functional is that applying if both the condensate and 
non-condensate fields were treated via a normally ordered characteristic func- 
tional. Here S(a;),S+(a;) are as in Eq. ()104p for the total field and 'ii{x),'ii^x) 
are the total field operators as in Eq.®, both involving all the modes and 
both of which are time independent. Details in deriving (jlOSp are set out in 
Appendix [9] Apart from the density operator the only quantity in the expres- 
sion (|107p for the characteristic functional that is time dependent is the factor 

exp 'l^^ J ^x Ec{x, t) E'^ [x, t)|-, and this simplifies the derivation of the func- 
tional Fokker-Planck equation. Eq. (|107l) relates the actual characteristic func- 
tional to that for the case where all modes are treated via a normally ordered 
characteristic functional. 

By comparing the expression ()107p for the characteristic functional with (|30p 
for the hybrid characteristic function it follows that the characteristic functionals 
and characteristic functions are equivalent. They both contain the same infor- 
mation - in one case as a functional of the fields Ec, S^t, Satc, S^^^ in the other 
as a function of all the expansion coefficients fo^' both condensate 

and non-condensate modes. 

x[Sc,S+,S^c,S+J=x(^,e+) (109) 

The phase space integral that relates the distribution function to the charac- 
teristic function can be expressed as a phase space functional integral in which 
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the distribution function P(a, , a* , a"*"*) is replaced by equivalent distribution 
functionals 

P[Tjj^,-ijy']=P{a,a+,a*,a+*) (110) 
where for short we introduce the notation 

The equivalence of functional and phase space integration based on separate 
modes can be established using the methods of functional calculus (see Appendix 
[HI see Ea. (l200p '). From Eq. pOl) the phase space integration is time independent 
even though the phase variables depend on time. Thus 

d^aciO) d^a-^iO) d^UNciO) d^a^^W (112) 

indicating that functional integration is also essentially non time dependent. 
Note that all modes are involved as are the four fields V'ci V'cj V'atC' V'tvc which 
are equivalent to the set of all ak , ■ Note that the functional integration 
J -D^V'c ^"^^ I ^^V'c involve a space grid with the same number of intervals as 
condensate modes. ). 

The characteristic function is still related to the quasi-distribution function 
P(a, a'^ , a* , a^*) via Eq. pip so writing this relationship in terms of functional 
integration we have an equivalent relation between the hybrid characteristic and 
distribution functionals 

xexp(i / dx{->pci^^i)'^ci^^'^)+-c{x,t)'ip'^{x,t))) 



xexp(i I dxV'jvc(2;'0'=ivc(^'*))-P[V^' V^*] ^'^P(* / dxENc{x,t) ^p%^{x,t)) 



X exp(i / dx {-ipci^^ "^(a;, t)+'ipffci^^ t) '^Nci^^ 0+ '^c{x, t) ^'^(a;, t)+ 'B.^cix, t) ipj^ci^^ 0)) 



X-PfV'o V'C' V'atO V'tVC: V'c. V'c^^WC) i^Nc] (113) 

where in the second line all the exponentials have been combined. The quasi- 
distribution functional is of the Wigner type for the condensate fields and the 
positive P type for the non-condensate fields, as in [14j . 
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These results are the same as those for the case of time independent phase 
space variables. We note that 



/ da;('0c(a;,i)S(t(a;,t)+V'7vc(2;,^)S^c(^'*)) = E "'^^fc 

keC,NC 
keC,NC 

which are time independent as before. Other details are given in Appendix 1^1 
3.3.1 Notation Change 

As in Section [5] we introduce the notation where we list the modes as 0^j,(x,i) 
with/i = -,+ and A = C,NC. Thus (t)Akix,t) =(I^Ak{x,t),(j)\^{x,t) =0^j,(a;,t). 
The condensate and non-condensate fields are listed as ip'^lxyt) and S^(a;,i), 

with ipci^' = "il^cix, t), tpci^^ = V'cCa;, t),i^Nci^^ = "^Ncix, t) and iptrci^^ t) 
ijj'jfQ{x,t). The total fields are 'ip^{x,t) = ijj^{x,t) + 'ip'^Q{x,t) and S^(x,i) = 
S^(a;, t) + E'^(-,{x, t). The couphng constants are listed as C^^. as in Eq. (|58|) . 

Writing {V'c' "^C' V'tvc V'^c) ^nd -0*.ee {i/-^, V'c*' V'atc' V'^ci t^en we 

have 

/// 1 D'^'^c D^i'c^'^^NC ^'^i^NC '^^ill be written as J D"^ ij; . In the new 
notation with phase space variables now a^Ak {t) and the field functions are 

tp'Xix, i) = E ^t^Ak{t)<l>Ak{x, t) (115) 

k 

Inverting Eg. dllSp gives 

"Mfe = / dx(j)^'^{x,t)^p'^{x,t) (116) 



Also from Eqs. ([55]), ([501) and (|116p we see that with 
d 



dyJ2'^Akiy^^^)(^AkBiKi('^''t) ^Aiy^*)+ / dyJ2^Aki'^^^)^AkBi'l'Bfiy^i) Kiy^t) 

kBl fcs; 
= - J d-vYl 'I^Bl(.^'''^)^BlAk(t>Akiy^t) + J ^yYl 'I^Ak(^'t)(^AkBl(t>Bliy^t) 



B^A B^A 



showing that the time derivative of the field function (C, NC) is a functional of 
both the original field function (C, NC) and the other field function {NC, C) . 
Note that the B — A terms cancel because of C^j.^; + C^ij^^. = 0. 



(117) 



39 



3.4 Quantum Correlation Functions 

The quantum correlation functions for the field operators can be obtained from 
the characteristic functionals via taking functional derivatives with respect to 
Sc, Sj, SjvCj "JVC ^^'•^ then setting these quantities to be zero. The result gives 
the quantum correlation functions as phase space functional integrals involving 
the distribution functionals. Thus the quantum correlation function is 

G(ri • • • Tp; • • • si; ; ui • • • • • • wi) 

= ({*c(?-l)^ ■ ■ • *c(7-p)^*c(Sq) • • • ^'c(si)}*Wc(wi)t • • • i'NcM'^'^NciVs) ■ ■ ■ *Arc(wi)) 

xV4(^i) ■ • • ^c{'rp)'4'c{sq) ■ ■ ■ V-clsi) V^jvcl^i) • ■ ■ V'^cK)V'jvc(«s) • ■ ■ i^NcM 

(118) 

For simplicity the t dependences of the condensate and non-condensate field 
operators and field functions have been left out. 

3.5 Functional Fokker-Planck Equations 

To derive the functional Fokker-Planck equation we first differentiate the expres- 
sion in Eq. (jl07l) for the characteristic functional with respect to t and second 
do the same to the expression (I113P involving the distribution functional. 
Differentiating the characteristic functional in its first form (|107p we have 



dt 



xTr(|expi j (*(a;) S+(a;))| (^^P^ jexpi j da;(S(j;) 



dt 



xr,.({»p.-/dx(*WH^W)} P {..,p./..(HW*.,x,)}) ,119) 

This result gives the time derivative of the characteristic functional as the sum 
of two terms. The first term in (jll9p only involves the time derivative of the 
density operator since from Eq. ()215p and ()42p the integrals involving and 
S+(a;) are time independent. This term equals the characteristic functional that 
would apply if p is replaced by -^p. The first term produces the standard terms 
in the functional Fokker-Planck equation. The second term depends on the 

time dependence of the / dx ^.(x, t) (x, t) factor due to the time dependent 
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modes and is equal to the characteristic functional multiplied by functionals of 
the condensate and non-condensate characteristic fields. The second term leads 
to additional diffusion terms in the functional Fokker-Planck equation. 



For the second term in M-x[^c^^'c^^NC^^tic\ quantity involving the 



time derivative is 

d r 1 



= 1^^^ y rf2^2c(a;,i)S(t;(x,t)| X exp|-i J dxEc{x,t)E^{x,t) 

(120) 

and using 

1 d f „ „ 
~2dt / ^^-c{x,t)^+{x,t) 

= U f fdxdy Iy.Y1 <l>Kx,t)Cki{t)My,t) \ {tENc{x,miE+{y,t)) \ 

['' [keCUNC J J 

+ i| / fdxdylY,J2 (bnx,t)CmMy,t) \ (*Sc(a;,t))(*S+^(y,t))l 

l"^ [leCkeNC ) ) 

(121) 

we find that the second term is given as the sum of two contributions as 

= f D^^P exp{i I dx'^^p^'{x,t)E-^'{x,t)) 

X / / dxdy \ \Y.Y. 
+ J D^i^ exp{i J dx ^^P''{x,t)E-f'{x,t)) 



§ s 

t)CckNClW'Pckiy>t) . ,+ I ^* 

JVC HNc\^,i)Hc\y^'^) 



I I dxdy U J2 Y^^iiMC, 

I keNC leC 



s s 

ciNCkit)(l>NCkiy^t)—^r—-———P[±,±*] 



(122) 



in terms of the new notation. In obtaining this result the exponential factor in 
Eq. (|120p is recombined via (|107p with the normally ordered characteristic func- 

tional Tri[ex,^|dxinx)EHx))} p {exp . / d.(.(.) ^t(..))|) to produce 

the original characteristic functional xi^c, Sj, ^nc, ^ncI- Writing the charac- 
teristic functional in the form in Eq. pi3p and noting that multiplication by 
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{iENc{x,t)){iE'X{y,t)) etc. can be replaced by functional differentiations such 
as ^ ^ — — result (|122p follows after applying functional integration by 

parts twice. Details for the derivation of (|12ip are set out in Appendix [HI 

For the first term in ^x[E!c, S^, S^vc", 5^^] we use the standard results 
for the correspondence rules associated with replacing the density operator 
by its product with the condensate and non-condensate field annihilation, cre- 
ation operators ilc{x,t),'^'^^{x,t),'^ Nc{x^t),ii'^^(j{x,t) that occur when is 
evaluated using the Liouville-von Neumann or master equation. For the non- 
condensate field operators we have 

P '^Nc{x,t)p x[Sc,S+,Sjvc,S+(g]^ X 

p => p'i'Nc{x,t) x[Sc,S+,SAfc,S+p] => ( .\ +iENc{x,t)] x 



5_ 

5{iENc{x,t))" 



For the condensate field operators we have 
p => '^cix,t)p X 
P ^ P^c{x,t) x[2c,5+,E;7> 
p ^ ^'^ci^-t)p x[Sc,S+,Sj. 



p =^ p^U^,t) 







1 , 
i. 


■Nc\ ^ 




2 




■i ' 


1 


wc] ^ 








■i ' 


1 


■Nc\ ^ 


\5{i3c{x,t)) 






■i ' 


1 , 
i, 


■Nc\ ^ 


\5{i^c{x,t)) 


2 



(124) 



From these results the corresponding changes to the distribution functionals can 
be deduced using (jlOSp and ()113p . This involves functional integration by parts. 
The correspondence rules for the non- condensate field operators are 

p ^ ^Nc{x,t)p P[^,-il^*]^i^Ncix,t)P 

P ^ pi'Nc{x,t) P[±^±*]^ (- , ,+ \ ,, +^Ncix^t)]P 

p ^]vc(^>i)P P[^,^1^f= \ +^lci^,t))p 
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The correspondence rules for the condensate field operators are: 



p ^ P^U^,t) PU,^*]^ (^+{x,t) + l—^—^P (126) 

where the i^c^'^c^'^C^i^c* ''^NC''^NC''^*NC^i^NC ^^'^ regarded as eight inde- 
pendent complex fields. The equivalence of these correspondence rules to those 
applying for separate modes can be established using the methods of functional 
calculus (see Appendix IH see Eq. ()188p which relates ordinary and functional 
differentiation) . 

By applying the correspondence rules (|125l) or (|126l) in succession to the var- 
ious products of the density operator with field annihilation, creation operators 
that arise from -^p we get for the first term in -^xi^Ci '^nC: -^jvc] 



J L>2^ exp(i J ^i/;^(a;,t)S-''(a;,t)) 



(127) 



in terms of the new notation, where A'^{x), D'^'^^{x,y) are the drift, diffusion 
terms associated with the standard functional Fokker-Planck equation. The 
integrals over x, y arise because the Hamiltonian involves spatial integrals over 
the field operators. Specific forms for the drift, diffusion terms are given for the 
cases of a single or two mode single component BEC in fl4l. 
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Differentiating the characteristic functional in its second form (|113p we have 



dt 



X exp(i / dx {-ipcix, t) t)+V'Arc(a;, t) '^Nci^^ 0+ '^c{x, t) t)+ Enc{x, t) ^pjjci^^ 0) 



d 

X Q^E'i^C: "00' '/'no V'wc V'C^C*' V'ArC-V'wc] (128) 



y exp(* J dxY,V{x,t)E-^{x,t))^^P[i^,^ 



in terms of the new notation. From Eq. piSp and (j42|) the integrals involving 
Ecixjt), SArc(x,i), and S^(-j(x,t) are time independent, whilst the 

functional integration J f f J D^ipc D^''l't;E''^'^NC ^^V'jvc b^ing the same as 

the phase space integration //// d^ac d'^a^d'^aMC d'^'^NC ~ II II ^^cic(O) d'^a^{0) d'^aNci^) <Pct.Nci^) 
- is also time independent. It follows that the time derivative of the character- 
istic functional is now determined from the time derivative of the distribution 
functional. 

Equating both sides of ^x[Sc, 2^, Satc, S^cl ^'i^- dH)) (|f27p and 
(|128p gives the functional Fokker-Planck equation for the hybrid distribution 
function. All the diffusion terms are combined, and there are new diffusion 
terms involving condensate fields paired with non-condensate fields arising due 
to the coupling coefficients. We have 



(129) 

where the new diffusion term is given by 

+ \ \ SAcSBNcSf,-t. i'^Yl^Aki^^t)'^AkBlit)'l>Bl{y,t)] \ 



k I 



+ i i^SscSANcS^-f, (^^'pBl{y^t)CBlAk(t)(t>Aki^^i)^ | 

(130) 

Details are given in Appendix |9l Here the additional diffusion contribution 
has been rewritten in a manifestly symmetric form. The drift A^(a;) vector 
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and the diffusion matrix E'^^g[x,y) are functions of the i/i^ and ultimately 
depend on x or x, y. The diffusion matrix is symmetric £'^'^(a;, y) = E'^'\{y, x). 
Thus for the hybrid case the functional Fokker-Planck equation has the same 
drift vector A'^{x) but a different diffusion matrix E'^'^g{x,y) to the diffusion 
matrix D^g{x,y) that would be obtained from the characteristic functional 
by applying the standard correspondence rules to the contribution. The 
additional terms involved the time dependent mode functions and the coupling 
coefficients. The functional Fokker-Planck equation can also be expressed in 
terms of real variables which involve the real and imaginary components of the 
phase variables tl)^, but we will not do that here. In deriving the functional 
Fokker-Planck equation there are often terms involving third and higher order 
derivatives arising from the contribution. These are usually discarded on 
the basis of being small due to scaling as higher powers of 1/^/N. 

The functional Fokker-Planck equation (jl29p for the distribution functional 
P[ , "0 *] must be equivalent to the ordinary Fokker-Planck equation Eq. (|63p 

for the distribution function P{a,a^ ,a* ,a^*). This is shown in Appendix [T] 

3.6 Langevin (Ito) Stochastic Field Equations 

To derive the Ito stochastic field equations we follow a similar approach to Sec- 
tion We consider a functional Flip] of the condensate and non-condensate 

fields and find conditions where the phase space functional average of the func- 
tional F[ip] is the same as the stochastic average when the fields ■0^4 are replaced 

by stochastic fields ip'^'' ■ 

3.6.1 Phase Space EXinctional Average 

The phase space functional average is given by 



For determining the quantum correlation function G{ri ■ ■ ■ Vp] Sq ■ ■ ■ si; ; mi • • • u^; Ws • ■ ■ ^i) 
the functional is = ^^(ri) • • • V'J (fp)0c(s«) ' ' ' ^ci^i) V'wc("i) ' ' ' 0Arc("'-)V'jvc(«s) ■ • • V'tvcC^i) 

- a product of functionals. In the phase space functional average both the dis- 
tribution functional and are time dependent, because from Eq. ()102|) the 

'i/'p, V'cj i^NC^ '4'~NC depend on time. However from (I112p the functional integra- 
tion is not time dependent. 

The time rate of change in the phase space functional average then consists 
of two contributions 




(131) 




(132) 
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To evaluate the first term we first consider changes d'ip'^{x,t) in tlj'^{x,t) 
during an interval St that results in a change (5F[V'] to the functional ^['0]. 

From Eq. (fT82| we find that 

so substituting from Eq. (|117p we get 

(134) 

Hence the first term in (jl32p is given by 



.2. d 



(135) 



Details are given in Appendix (TU] 

Using the functional Fokker-Planck equation (|129p we find after using inte- 
gration by parts that the second term in (|132p is 

(136) 



where functional integration by parts has been used. Details are given in Ap- 
pendix [TOl 
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Overall we now have 





JEE / f'^^'^y 



fj.A yB 



This result can be compared to that for a function f (a, a"'") as set out in 
Section [2] where from Eq. ((70|) 



Si 




F{a,a+) 



^Ak Bl 



(138) 



We see that the derivative of the functional phase space average is not the same 
form as that for the original phase space average involving separate modes, 
since a term corresponding to one of the first order functional derivative terms 
is missing. As a specific case consider the functional F[^] = iljQ{ri,t)^(j{si,t) 

that is involved in working out the symmetrically ordered correlation function 



(137) 



({*c(ri,t)t^'c(si,i)}). By substituting we get = 



— )• 



E^CfcW-^Cfc 

k 



I 



hi 



'Pcki^ii^) '^ckW'^ciit) 0ci(*i'^)- The last equation is useful for work- 



ing out the time derivatives of the various F{a,a'^) — a^jj^aci, but it does 
not include the contribution to the time derivative of F[^] from the various 

</'cfc(^ii ^)<^Ci('*i' These contributions can of course be worked out and in- 
cluded later, but clearly the advantage of using field functions rather than sep- 
arate phase variables is in terms of focusing directly on the quantum correlation 
functions that involve condensate and non-condensate field operators rather 
than just mode operators. As we will see there is no paradox associated with 
the missing term - there is a compensating term arising from the stochastic field 
average. 
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3.6.2 Stochastic Fields Average 

For the Ito stochastic field approach we now replace the fields tp'^ by stochastic 
fields iPa"^- We suppose ipA^ix^t) satisfies an Ito stochastic field equation of the 
form 

/t+st 
dtiTaiti) (139) 
u 

or equivalently 

^^r/{x,t) =g^jx) + Y,M^^ix)r,it+) (140) 

a 

where expansions of tp'^^ {x,t),g'^{x) and Af^a (^) terms of the mode functions 
(j)'^j.{x,t) and stochastic phases ^^^^(i) will be determined later. 

Using an obvious generalisation of (I189P we can write the change in the 
stochastic functional Flip^ , tp''^^] due to changes (5'0^*(a;, t) in the stochastic fields 
as 



F[^'{x, t) + S^j'ix, t),ij'+{x, t) + 6ij'+{x, t)] - F[^'{x, t), ^"+(x, t)] 



+ IJJ dxdy SrA-, t) Srs'iy, t) ^ ^ 

correct to second order. Using the result that for any function H{ip^ [x^t)) of 
the stochastic fields we have 



H{,l;^'j^{x,ti))^a{t2)Tb{t:i)V,{t^)..Tk{ti) 



= Hiij^'^ixM)) ^a{t2)n{t3)T,{ti)..Vk{tl) ti <t2,t3,..,ti (142) 

we can then evalua,te Flip' {x,t) + dip" {x,t),il}'+ {x,t) + Sip"^ {x,t)] - F[ip" {x,t),ip"+ {x,t)]. 

For the first order derivative terms we have after substituting from (|139p . 
expanding and carrying out the stochastic averaging 



/ dx \j2srAix,t) 



(SFlip^iP 



V SlpAi^^'t) 



- /-E(M)/5<^)« 

where the stochastic average rules for sums and products have been used, the 
non-correlation between the averages of functions of ip"^(x^t") at time t and the 

r at later times between t\,ot-\-bt\% applied and the term involving VaitLx) is 
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equal to zero from (1751) . Note that this term is proportional to St. Details are 
set out in Appendix [TOl 

For the second order derivative terms we have after substituting from (|139p , 
expanding and carrying out the stochastic averaging 



i 1 1 dxdy J2 srA{^,t)srBiy,t) 



dxdy 



s+l 



St 



dxdy 



liA.uB 



s+l 



Sil}'X{x,t)Sil}''B{y,t) 



Wix) M^{y)rA. 



A,B 



St 



(144) 



where the stochastic average rules for sums and products have been used, the 
non-correlation between the averages of functions of ?/'^(x,t) at time t and the 
r at later times between t to t + St is applied. The terms involving a single T 
have a zero stochastic average, whilst from Eg. ([75]) the terms with two F give a 
stochastic average proportional to St. In evaluating the latter term Eq. ((8T|) is 
used. Details are set out in Appendix ITOl 

The remaining terms give stochastic averages correct to order St^ or higher 
so that we have correct to first order in St 



F[ib'{x,t) + Sij'ix,t),i/+{x,t) + S^'+{x,t)] - F[^%x,t),iP'+{x,t)] 

St 



fdxY^ 



fSF[iP',i(;' 



s+l 



\ Sil}'^j^{x,t) 



(145) 



or 



1 



d_ 

dt 



dxdy 



hAmB 



s+l 



Ha ix,t)SiJyiy,t) 



Mx) A/-^(2/)]rB 



St 



F[rix,t),r^ix,t) 



[dxY: 

•' fiA 



fSF[%p'',tl)' 



s+l 



(146) 



1 



dxdy 



5'F[i}',il: 



s+l 



[Af{x) U^{y)]rB 



For the phase space functional average (1137^ and the stochastic average (I146P 
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to agree we require two conditions to be satisfied, as follows. 




W{x)M^{y)]XB^E'X's{x,y) 




(147) 



k I 



+ 2 \ SbcSancS^-^, i'^'^(t)Bl{y,t)C^iAk{t)4>Aki^^t) 



k I 



These conditions determine the quantities G^Ai^) ^^'^ ^Aai^) that occur in the 
Ito stochastic field equations (I139P in terms of the quantities that occur in 
the functional Fokker-Planck equation (|129p . Note that the B = A terms in 
the result for ^^(a;) actually cancel out because C^^.^; + CAi^Ak ~ ^- These 
relationships are to be interpreted as replacing the stochastic quantities tj^(a;), 
■^Aai^) the Ito stochastic field equation (|140p with the quantities A'^{x), 
Nj^^{x) and tl)'^{x^t) originating from the functional Fokker-Planck equation 
(|129p . which are then regarded as depending on the stochastic fields ip^^{x,t). 

We next show that these conditions can always be satisfied. In Appendix 
[7] we have shown that the Fokker-Planck equation for the distribution function 
and the functional Fokker-Planck equation for the distribution functional are 
equivalent, and we can make use of the results obtained there. 

First, we consider the condition for Q''^{x). From (|203p the drift term A'^{x) 
in the functional Fokker-Planck equation is related to the drift vector A^j, in 
the Fokker-Planck equation via 




Ak 



k 



^Ak 




(148) 
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We can also use Eqs. (|115p . (|116p and (pO]) to show that 

Bl't'Bl 

iy,t) 

B \ kl / 

B \ kl / 

= E'^Afc(2;,0E'^Afci3i"MB' + E 
\ fc Bl J \ k 



^^rAki^,t)a^Ak] (149) 



Clearly the first condition is satisfied via the choice 



g^{x) ^ E'^Afc(^'*Kfc + E'^Afc(^,OEC^fcB/«;.B/+E|'^Afc(^'^)"; 



k k Bl k 

(150) 

written in terms of mode functions and quantities derived from the Fokker- 
Planck equation 

Second, for the condition for M'^g^{x) we need to be able to find M'^^{x) such 
that 

E-^Aa(^)A/'Ba(y) ^ E^^^Bi^^y) (151) 
a 

Using the relationship in Eq. pOSp between the diffusion matrix Ej^'^g{x,y) in 
the functional Fokker-Planck equation (|129l) and the diffusion matrix E"^^ in 
the Fokker-Planck equation (|63p we have 



EAsi^^y) = Y.'f'Aki'^^^^^AkBl^t^Bliy^t) 
kl 

^AkBi ^ dxdy (p^'i^{x,t)E'^''s{x,y)4>B'j^ {y,t) (152) 



The matrix E is symmetric {Ej^i^^i — -E^'^^^,), so from the Takagi factorisation 
[45] . [46] it can be written as KK^, where the matrix K is as in Section [2j 
Eg. (|86l) . Thus the diffusion matrix E'^'^g{x,y) may be factorised as 

ETB{^.y)^Y.^Aai^)^Ba{y) (153) 
a 

where 

Ka{^) = Y^^Aka^AkM (154) 

k 

ETkBi - T^^AkaKma (155) 



51 



Hence the second condition is satisfied with the choice 




(156) 



fc 



again written in terms of mode functions and quantities derived from the Fokker- 
Planck equation 

Finally, we can show that the Ito stochastic field equations (|140[) for the 
stochastic fields ^p'^ (x, t) in the field theory treatment are entirely equivalent to 
the Ito stochastic equations (I57|) for the stochastic phase variables a^^^, in the 
separate mode treatment. The proof is given in Appendix [T] 

This proof shows that the stochastic fields are given by 



where the stochastic feature results only from the phase variables being replaced 
by stochastic variables. Thus the stochastic field equations for i/jj^{x,t) are 
equivalent to the stochastic phase variable equations for the 

3.7 Expressions for Ito Stochastic Field Equations 

The original Ito stochastic field equation (|140p can also be written in terms of 
the drift, diffusion terms that are derived directly from the functional Fokker- 
Planck equation (|129p itself. Substituting from (|156p and (I147P we have for the 
Ito stochastic field equation 




(157) 



k 






(158) 



a 
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where Aj^^x) is the drift term and N'^^{x) is related to the diffusion term 
Ej^"g{x,y) that occur in the functional Fokker-Planck equation (I129P via 

a 

+ 2 \ ^AC^BNC 



^.-^ EE 

t)C^AkBlit)^Bliy,t) 

\ k I 

+ n i ^BC^ANC Su^u [T.J2^Bliy,t)C'^,lAkitWAki^,t) 



The B = A terms have been cancelled out using C^j,^; ^^AkAi ~ ^- "^^^ result 
for ■^ip'^^x, t) may be compared to the Ito stochastic equation ((87|) for the phase 
variables. The term that is added to yl^(x) involves a non-local contribution 
with a kernel that depends on the mode functions and their time dependence via 
the coupling coefhcients. It has the effect of coupling the condensate stochastic 
fields to the non-condensate stochastic fields V'tvc = ~7 +)• This term is 
not present when time independent modes are used. Note that it is of the same 
form as in (|117p for the time derivative of the (non-stochastic) field functions 
'(/'^(Xjt). This form for the Ito stochastic field equation is the most useful, since 
it involves the drift ^^4(2;) and diffusion terms D'^'^g{x, y) in the standard func- 
tional Fokker-Planck equation obtained by applying the correspondence rules for 
the contribution in the characteristic functional, plus the time dependent 
modes (/)^j,(a;,t) and the coupling coefficients C^^gi{t). 

3.8 Classical and Noise Field Terms 

We can write the Ito stochastic field equation in terms of a classical term and 
a noise term 

^i^aM = (^^^rA{^,t))^^^^+{^l-^r/{^,t) 



' class 



- [dyJ2 (T.^Aki^^^)(^lkBi'I^B^iy^t)] rB{y,t) 

B^A \ kl ) 

J dy E (E<^Al^(y'^)^Afc%'/'^;(^>o) VAiy.t) 



B^A \ kl 



^^rAi^,t)) = Y.^Aai^)^aiU) (160) 



noise 



If only the classical terms were included, then the solution for tp'^{x,t) would 
determine classical field functions; hence their name. Note that the number of 
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noise terms will be equal to 2{nc+nNc), or twice the total number of condensate 
and non-condensate modes. 



3.9 Properties of Noise Fields 

The stochastic averages of the noise field terms can now be evaluated. For a 
single noise term and the product of two noise terms we have using Eqs. (|78p 
and (fT42ll 



showing that the stochastic average of a single noise field term is zero whilst 
that for the product of two noise field terms is delta correlated in time and equal 
to the appropriate diffusion matrix element. Note that it is not necessarily delta 
correlate in space. 

In fact the diffusion matrix elements determine all the stochastic averages 
of products of noise field terms. With an odd number of terms the stochastic 
average is zero. For an even number of terms the stochastic average involves 
sums of stochastic averages of products of diffusion matrix elements, which 
reflects the Gaussian-Markoff properties of the Ta- Thus for three and four 
noise field terms 

^^m(^I'^i)) f^V'^s(a;2,t2)j f^V'5c(a^3,i3)j =0 (162) 



^i'U(xi,ti)^ (^^CB(a;2,i2)^ (^■^i^lcix3,t3)^ (^^V' AD (2^4,^4)^ 
= S{ti ~ t2)6{t3 - ti) E'^g{xi,X2;ti,2)E^')j{x3,X4,;t3^4) 



+6{ti - t3)S{t2 - h) E%{xi,x:i-ti^3)El^^{x2,Xi-t2A) 
+6{ti - ti)6{t2 - h) E'^^^{xi, X4;tiA)E''J^{x2, X3;t2,3) (163) 
Note that the noise field terms are not themselves Gaussian-Markoff processes. 
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4 Conclusions 



This paper has sot out a phase space theory approach to treating the dynami- 
cal behaviour of Bosc-Einstein condensates which could be applied to a variety 
of experimental situations, inchicling interferometry with BEC involving time- 
dependent double well potentials. The key feature is the use of time-dependent 
mode fimctions, chosen so that a few (one, two,..) highly occupied modes pro- 
vide a good approximation to the physics of interacting bosons in the time 
dependent potential at temperatures well below the BEC transition temem- 
pcrature. The motivation is that when the effect the lightly occiipiecl non- 
condensate modes are taken into account in treatments of decoherence effects 
allowing for Bogoliubov excitations, collisions with thermal bosons etc., the 
previously determined condensate modes are still important in treating these 
extensions. To allow for time dependent mode functions, the major innovation 
here is the introduction of time dependent phase variables, the time depen- 
dence being chosen to match that of the time dependent mode annihilation, 
creation operators. This means that in this approach we are representing time 
dependent mode annihilation, creation operators by time dependent phase vari- 
ables, whilst time independent total field annihilation, creation operators are 
represented by time independent field functions. In the more usual approach 
using time independent phase variables the opposite situation applies, which 
is rather unsatisfactory. The theory presented here treats the two situations, 
one {mode theory) being where mode annihilation, creation operators and their 
related phase variables and distribution functions are dealt with specifically, 
the other {field theory) being where field creation, annihilation operators and 
their related field functions and distribution functional involve a description 
where individual modes are not distinguished. Though each situation is treated 
separately, they are shown to be equivalent. Within each of these two situ- 
ations the modes are divided up between condensate (highly occupied) modes 
and non-condensate (sparsely occupied) modes. This is referred to as the hybrid 
approach. For the hybrid case a Wigner type distribution function or functional 
is used for the condensate bosons and a positive P type for the non-condensate 
bosons. Fokkcr-Planck or functional Fokker-Planck equations are derived, along 
with the Ito stochastic equations for stochastic phases or stochastic fields, and 
the relationship between the Ito equation quantities and those in the Fokker- 
Planck equation is obtained. The stochastic equations contain the sum of a 
classical term related to the drift vector and a noise term, and the stochastic 
properties of the noise terms are related to the diffusion matrix. Also, the ex- 
pressions for the stochastic fields involve the same expansion in terms of time 
dependent mode functions and stochastic phase variables, as had applied to the 
field function expansions involving non-stochastic phase variables. 

The key results are set out below. For the hybrid approach the condensate 
and non-condensate field operators are time dependent, as are the correspond- 
ing field functions. Also, both the Fokker-Planck and functional Fokker-Planck 
equations differ from those derived using the usual correspondence rules, the 
drift vectors are unchanged but the diffusion matrices contain additional terms 
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relating to the time dependence of the mode functions and involving the coupling 
coefficients, which are defined via overlap integrals between mode functions and 
their time derivatives. In addition, there are extra terms in the Ito stochastic 
equations both for the stochastic phases and stochastic fields, again related to 
the time dependence of the mode functions and involving the coupling coeffi- 
cients. For the hybrid case the final form for the Ito stochastic field equation 
in Eqs. (|164|) . (|165|) is the most useful, since it involves the drift A'^{x) and 
diffusion terms D^g{x,y) in the standard functional Fokker-Planck equation 
obtained by applying the correspondence rules for the contribution in the 
characteristic functional, plus the time dependent modes and the cou- 

pling coefficients C^^ {t) . 

From the point of view of best describing the physics of BECs even near zero 
temperature the treatment preferred here is the hybrid one involving condensate 
and non-condensate fields rather than considering a large number of separate 
modes, but based on time dependent mode functions. For this reason, the field 
theory treatment has been given an important emphasis in this paper. That 
treatment of the hybrid approach capitalises on the simple but accurate picture 
based on highly occupied condensate modes and lightly occupied non-condensate 
modes, but avoids a mode-by-mode approach in determining the quantum cor- 
relation functions - which after all only involve field operators. There are of 
course some costs, for numerical work the set of time dependent mode functions 
and the coupling constants would need to be worked out first, but hopefully 
only a few condensate modes would be needed and the non-condensate modes 
might be introduced via physically based variational methods involving Fock 
states where only one boson is in a non-condensate mode, or by just applying 
Schmidt orthogonalisation procedures. 

In Ref. [2] the hybrid approach was applied based on including time de- 
pendent condensate modes and the functional Fokker-Planck equations derived 
using the correspondence rules as set out here in Eqs. (jl25l) and (|126p . In terms 
of this paper (see Eq. (|166p ) this gives the drift vector and the diffusion 

matrix contribution D'^'^g{x,y). However, no terms were included in either the 
functional Fokker-Planck equation or the Ito stochastic field equations involving 
the coupling coefficients C^^, (t) . The theory in Ref. [13] was based on using 
time independent phase variables and characteristic function variables, but the 
time dependences of both the condensate and non-condensate field operators 
and field functions consequential on using time dependent condensate modes 
were ignored. Eqs. (192), (193) and (207) in Ref. [14] do not include contribu- 
tions involving the coupling coefficients. The results for the Ito stochastic field 
equations in Ref. [Ij are correct if time independent (or slowly varying) con- 
densate and non-condensate modes were used, but then the condensate would 
require many more than one or two modes for an adequate description. Note 
however that the non-local term involving F{r, s) in the diffusion matrices in 
Ref. [13] would still be present, as this arose from the -^p contribution in the 
characteristic functional and does not involve mode time dependence. However, 
results from the present paper are really needed to supplement those in Ref. [T4] 
if the physically based idea of just having one or two genuinely time-dependent 
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condensate modes is to be implemented. 

The results in this paper can be easily converted back to the corresponding 
results based on using a convenient set of time independent mode functions. In 
this case all the coupling constants become zero and the phase variables, mode 
and field operators become time independent again. The functional Fokker- 
Planck and Fokker-Planck equations would revert to those obtainable by stan- 
dard methods using the correspondence rules, the Ito stochastic equations for 
fields and phase variables are also the same as in standard methods, as also 
are the relationships between these equations and the Fokker-Planck equations. 
However, such mode functions are not closely aligned to those that provide a 
good approximation to the near zero temperature physics, so any advantages in 
terms of simplicity in applying the theory and convenience for numerical work 
is likely to be negated. 

As well as the hybrid approach treated in this paper, situations where time 
dependent mode functions are involved may also be treated using the so-called 
combined approach, where the modes are not divided up into condensate or 
non-condensate types and only the total field operators and field functions are 
considered. For completeness, the results for the combined approach are also set 
out below. The same method involving time dependent phase space variables 
can be used. Again, two situations may be considered, one (mode theory) being 
where mode annihilation, creation operators and their related phase variables 
and distribution functions are dealt with specifically, the other (field theory) 
being where total field creation, annihilation operators and their related field 
functions and distribution functionals involve a description where individual 
modes are not distinguished. The distribution function or functional may be 
either of the Wigner or positive P type. For the combined approach both the 
Fokker-Planck and functional Fokker-Planck equations are exactly the same as 
those derived using the usual correspondence rules, even though time dependent 
modes are present. The field functions are time independent. However, there are 
extra terms in the Ito stochastic equations for the stochastic phases though not 
for the stochastic fields. These terms relate to the time dependence of the mode 
functions and involve the coupling coefficients. For the combined case, the final 
form for the Ito stochastic field equation in Eqs. (|17ip is the most useful, since it 
involves the drift A^^ (x) and diffusion terms D^^ '^{x,y) in the standard functional 
Fokker-Planck equation obtained by applying the correspondence rules for the 
contribution in the characteristic functional, plus the time dependent modes 
(j)'^{x,t) and the coupling coefficients C'^i{t). 
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5 Summary of Key Results 

5.1 Hybrid Approach - Condensate and Non- Condensate 
Modes and Fields 

Here we summarise the key results obtained in this paper. Resuhs for func- 
tional Fokker-Planck, Ito stochastic field equations and the Fokker-Planck, Ito 
stochastic phase variable equations for the hybrid approach where the time de- 
pendent modes are divided into nc condensate modes and njvc non-condensate 
modes, plus conversions and key results are as follows. 
For the field theory results: 



k ■' 

^^r/{^,t) = A'X{x)+ jdyY, (Y.^Ak{^^^)c'XkBi<t>B't{y,t)^ rB\y,t) 

+ IdyY. \ T.^~A^MC-A^Bi4>''Bi{^,t) \ rA{y,t) 

•' B^A \ kl J 

+ E^l(^)r«(*+) (164) 
E^^''^ix,y) = D^/^{x,y) 

+ ^ jc^AC^B Arc (5^-1. (^J2'^Aki^^t)C'^,^ Si{t)(j)Bi{y,t)^ | 

+ ^ ^^SBcSANcSu-f, (^Y^(l)Bl{y,t)CBi Ak{t)(I^Ak(^^t)^ | 

= Y^Kai^)Kaiy) (165) 

a 

rAix,t) = ^a^^AkmikM c,;^Ak= [ dyci>Tk{y,t)rA{yMm 



58 



For the mode theory results: 



d 



= T.^AkaK'hla (167) 

a 

Conversions between the quantities are as follows 

k 

^AkBi = / dxdy(l)^'^{x,t)E'XB{x,y)(f>B'^{y,t) E^X^ix^y) ^^<l)^Xki^,t)E'^l gi(j)Bi{y,t) 

kl 

Kai^) = Y.^Aka<^AkM (168) 

k 

Key results for the time dependences are 
d 

Bl 

d 

;af,Ak = "^C^kBiO^t^Bi (169) 



dt 

with coupling coefficients 



Bl 



CAkBlit) - j dx^^^^C^si{x,t)^CAkBl[t) 

ClkBiit) = ldx^^^^cj^si{^,t)^C*A,Biit) (170) 

5.2 Combined Approach - Modes and Fields 

A much simpler treatment applies if all n condensate and non-condensate modes 
are considered together in the so-called combined approach, even though the 
mode functions are time dependent. Here we will just summarise the results 
for completeness, but without proof since the proofs are easily obtained from 
the previous sections. Many of the extra terms for the hybrid case just become 
zero. The phase variables are again considered to be time dependent, as in Eq. 
([37)) as are the characteristic function variables, as in Eq. (|411) . The phase space 
treatment is based on the positive P distribution for all modes, but similar re- 
sults would apply if the double space Wigner representation was used. It then 
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turns out that the Fokker-Planck equations involving the phase variables a^k 
{fi = —,+,k = 1, ..,n) have the same drift A'^ and diffusion D'^^ terms as in 
the standard Fokker-Planck equations resulting from just the contribution 
in the characteristic function. However, the Ito stochastic equations have an 
additional contribution J^i '^fii to the classical field term, though the noise 
term just involves the matrix B related to the diffusion matrix via D = BB^ . 
The functional Fokker-Planck for the time independent total fields ip^ {x) turns 
out to be exactly the same as in the standard functional Fokker-Planck equa- 
tions resulting from just the contribution in the characteristic functional, 
involving drift A^{x) and diffusion D'-"^{x,y) terms. Also, the Ito stochastic 
field equations are of the same form as if the mode functions were time indepen- 
dent. The classical field term is given by the drift term and the noise field term 
is related in the standard way to the diffusion term in the functional Fokker- 
Planck equation. Furthermore, the stochastic field tp''^{x,t) can be expanded in 
terms of the time dependent mode functions with the stochastic phase variables 
as coefficients, with the same form as for the non-stochastic fields 'ip^(x). 

Results for functional Fokker-Planck, Ito stochastic field equations and the 
Fokker-Planck, Ito stochastic phase variable equations for the combined ap- 
proach where all time dependent modes are considered together, plus conver- 
sions and key results are as follows. 

For the field theory results: 





k 





(171) 



a 



a 




k 




(172) 
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D'i;{)P{a,a+,a',a^*) 



For the mode theory results: 

/ a 

Dti = (173) 

a 

Conversions between the quantities are as follows 

^ k 
Dti = j [dxdyct>^^{x,t)D^-'{x,y)cl>r{y,t) D>^%x,y) = Y,€{^,t)Dti 4>'{{y,t) 

kl 

N^{x) = ^Bir^M (174) 

k 

Key results for the time dependences are 
d 

df 

with coupling coefficients 

cm = j dx^-^^cl>*i{x,t)=Cm (176) 



y^ki^^t) = E^«'<^r(^>i) 

I 

d 
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6 Appendix - Functional Calculus 



6.1 Functionals 

Essentially a functional F[ip{x)] of a field function tp{x) just defines a process 
that results in a c-numbcr which depends on all the values of the field func- 
tion, that is over the entire range of positions x. Unless otherwise stated the 
field function is complex. Generalisations to cases involving several functions 
'tp{x),^~^ (x) or il;(x),^~^ [x),ip* {x),ip~^* (x) etc arc straightforward, as also when 
X refers to position r in a 3D space. The field functions may also sometimes be 
time-dependent. 

Examples of functionals are numerous. If x(a;) is fixed then the integral 

J dxx*{x)ip{x) is a functional of ip{x), often written xbPi^)]- This functional 

defines the scalar product, often written (x, V^). Herex*(x) acts as a kernel. A 
function 'ip{y) can also be considered as a functional of il^{x), with the delta 
function acting as the kernel. This specific functional is written as Fy[il){x)] = 

j dx5{y - x)ii){x) = tp{y). 

If the function il>{x) is expanded in terms of orthonormal mode functions 
(j>k{x) with expansion coefficients 

i,{x) = Y.'^Mx) (177) 

fc 

c^k = j dx<l)l{x)ij{x) (178) 

then the functional F[iIj{x)] can always be determined from the ak, so is the 

same as the function /(ai, 0:2, ••, ctk, ••) obtained by replacing ip{x) by its mode 
expansion in the process that determines the functional. 

F[V(x)] =/(ai,a2,..,afc,..) (179) 

The orthonormality and completeness relationships are 

dx<j)lix,t)<j)i{x,t) = Ski (180) 
Y,M^,t)rk{y,t) = S{x-y) (181) 

k 

6.2 Functional Differentiation 

A summary oi functional differentiation is as follows. For a functional Fltp^x)] 
of a field the functional derivative ^"^2)^ ^ defined by 

F[i,{x) + Sij{x)] = F[i,{x)] +Jdx 5^{x) (^^^^) (182) 
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where dip{x) is any small but arbitrary change in ip{x). In this equation the 
left side is a functional of tl^{x) + 5ip{x) and the first term on the right side 
is a functional of ipi^)- The second term on the right side is a functional of 
Sipix) and thus the functional derivative must be a function of x; hence the 
subscript x. In most situations this subscript wih be left understood. If we 
write Sip{x) = eS{x — y) for small e then an equivalent result for the functional 
derivative at a; = y is 

f6F[i^{x)]\ _ ^.^^^ / Fji^jx) + eSjx - y)] - F[i;{x)] \ ^ ^^^^^ 



V J x=y 

Note that for a fixed function ^{x) 

sax) 



(184) 



S'>p{x) 

since the fixed function does not change when ^'{x) does. On the other hand 

Functional differentiation satisfies many of the rules of ordinary differentia- 
tion including a product rule 



Note that for a fixed function ^x) we have from Eqs. (|186p and (|184p 

5{ax)Fmx)]} f SF[r^ix)] \ 

6^{x) -^^""'y Si;{x) J ^ ' 

so a fixed function may be moved through the functional differentiation. 

Mode expansions can be used to relate functional differentiation and ordinary 
differentiation. 

f 5F['4;{x)]\ ,*/ . 9/(ai,a2, ..,afc, ..) 



k 



dak J \ Stpix 

Note that the first of these simple results are obtained by considering in (|182l) 
the change in Sip{x) ~ ''^^Sak 4'kix) given by changes in the ak, and then 

fe 

equating F[ip{x)+S'ip{x)] — F[ip{x)] to the corresponding result f{ai+dai,a2 + 
da2,--,ak + dak,--) — /(ai, 02, .., afc, ..). The second is then obtained via the 
orthogonality result (I180p . 
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A Taylor series expansion to higher orders can also be obtained. Thus to 
second order 

(189) 

and can be obtained via a Taylor expansion of /(ai,a2, ■■,ctk, •■) plus the rela- 
tionships in (|188p . 



6.3 Functional Integration 

A brief summary of functional integration is as follows. If there are n modes 
then the range for each function ip{x) is divided up into n small intervals Axi = 
Xi+i — Xi (the iih interval, where e > |Axi|), then we may specify the value ip^ 
of the function ip{x) in the ith interval via the average 

^^ = ^ I dx^ix) (190) 



A, 



Xi 



and then any functional F[ip{x)] may be regarded as a function F(iIji,iI>2, --jipi, --jipn) 
of all the t|;^, and ordinary integration over the -0^ is used to define the func- 
tional integral. If each function ip{x) = "ipxi^) + iipy{x)As written in terms of its 
real and imaginary parts, then the functional integration becomes an ordinary 
integration over the values tp^^ , ifj^y of these components in each interval i of the 
function F{ipi,ip2^ ■■y'4'n) multiplied by a suitably chosen weight function 
iy(?A]^, ?A2: - -I "011 V'n)- Thus the functional integral is defined by 



D ■>jjF['ip{x)] = lim lini / ••• / d ■02--c^ "^j-c^ V'„'«('0i,V'2>-->i,--,V'„) 

xF(V'i,V'2,..,V„..,V'„) (191) 

where the number of modes is increased to infinity along with the space interval 
decreasing to zero. The symbol D^tp stands for (Pipi(P%p2--d^^i--d^^n ^(V'l: V'2j V'i: --j V'n)' 
where the quantity d^V'i means dip^^dipiy 

A useful integration by parts rule can often be established from Ea. (|186p . 
Consider the functional H[il}{x)] — F[^{x)]G[il;{x)]. Then 

(192) 

Functional integration and phase space integration are inter-related via the 
mode expansion, so we can relate the value ip^ of the function in the ith interval 
to that of the mode function (pi.^^ via 

tA, =^afc0fc, (193) 

k 
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Choosing the number of intervals to coincide with the number of modes the inte- 
gration J ' ' J d^V'i'^^V'2--'^^^i--'^^''/'n '^^'i tie changed to an integration J ' ' J (Poti'Pc(2--(Poik--(P 
over the expansion coefficients ai, 02, .., ak, --an- We have 



D^iJF[^P{x)] 



hm hm 

n— foo €— 



(fai(fa2..(fak--d'^an || J(q!i, 02, .., afe, ..a„)|| 
xv{ai,a2, ■■,ak, ..a„) /(ai, 02, ..,afe, ..a„) (194) 



where w(ai, 02, .., ak, --an) is the function that equals w{ipi, ip2, ■■, V'i: --j V"™) and 
the Jacobian is given by 



II J(ai, a2, .., ak, ■ 



Now using Ea. (fT93)) 



dakx 
dakx 









9^1^ 


fV-lx 


av-ix 






ga2x 








da2y 


dO-ny 








3V'2x 


3V'2x 


9V'2x 


'9'02x 




9a2x 




SOnx 




da2y 










9'/'„x 


fi'V'^.x 




av-^x 




9a2x 




9a7ix 




da2y 














di'ly 




gaix 






aa„x 


daiy 


dOL2y 


da^y 




d4>2y 






d4'2y 


9i'2y 


di^2y 


9aix 






flanx 


daiy 


da.2y 


9a„„ 




9^'^y 




d'^'^y 




d^^y 




9aix 






Sq„x 


daiy 


da2y 
















(195) 
















= 




Oaky 


= -0. 


iy 







l^kiy 



^^^y 



dak 



ky 



(196) 



and evaluating the Jacobian after showing that {J J'^)iakb = ^ik^ab/ using 
the completeness relationship in Eg. dlSip we find that 



||J(ai,a2,..,afe,..a„)|| = IT (Ax^ 



(197) 



and thus 



lim lim 

n— 7-00 e— ?-0 



dPai<Pa2--d?ak..d^an TT , . ^ , 



xv{ai,a2,..,ak, •■«„) /(ai, 02, .., afc, ••«„) 



(198) 



This key result expresses the original functional integral as a phase space integral 
over the expansion coefhcients ak of the function ^{x) in terms of the mode 
functions (pkix). 

The general result can be simplified with a special choice of the weight 
function 

iz;(V'i,V2,-,V'„-,V'„) = n(^^^) (199) 
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and we then get a simple expression for the functional integral 



/ D'^ip F[ip(x)] = lim lim / ••• / (fai(fa2..(fak--(fanf(ai,a2,..,ak,--an) 

J n^ooe^Oj J 

(200) 

In this form of the functional integral the original functional F[iIj{x)] has been 
replaced by the equivalent function /(ai, a2, ••, Q^fc, ••Q:„) of the expansion co- 
efficients ak, and the functional integration is now replaced by a phase space 
integration over the expansion coefficients. 
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7 Appendix - Equivalence of Separate Mode and 
Field Theory Treatments 

7.1 Fokker-Planck and Functional Fokker-Planck Equa- 
tions 

The functional Fokker-Planck equation for the distribution functional P[V-' , "0 *] must 

be equivalent to the ordinary Fokker-Planck equation Eq. for the distri- 

bution function P(a, «+, a*, a"*"*). The drift and diffusion terms can be inter- 
related because functional and ordinary differentiation are related, as explained 
in Appendix ini We use 

d 



^ dxc^'X,{x,t)(-J—^\ (201) 



in the functional Fokker-Planck equation (I129P which becomes 

= |-P(a,a+, «*,«+*) 
ot 



I Ai I 

y fiA k 



a , a ) 



' E ^^A.- + I E E ^ ^9;^^^^^^^' 1 «^*) (202) 



and gives the Fokker-Planck equation (1551) with 

A'Xi^) = E'^A.(^,^)^i^fe 



^AfcBi = y y dxdy(P^'i^{x,t)E'^''s{x,y)(f)s'^{y,t) 
ETsi^^y) = Y.'t^AkME'AkBi^miy.t) (203) 
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which is the same as Eqs. ([M]). Note that in particular 

-"Am Bn 



J J dxdy4>^'^{x,t) ^^^^Ac6BNcS^,-u (^J2^Ami^^t)(^'A 

+ J J dxdy<j);^'^{x,t) ^hscSANcS^^f, ij2J2'^Bniy^*)^BnAmit)<pAmix^t) 



(t'B^iiy^t) 



{^AC^B NC Sfj.-^C'^f^ + Sbc^ANC K-fiCsij^i,} (204) 



_ 1 

2 

as in Eq. (IMl) . Thus the functional and ordinary Fokker-Planck equations for 
the hybrid case are consistent with each other. 



7.2 Ito Stochastic Equations for Fields and Phase Vari- 
ables 

To see how the stochastic fields i/jj^{x,t) are related to the stochastic phases 
^^Afc' ^^'^ that the Ito stochastic field equation (|139p becomes on substitut- 
ing for gf'ix) and Af^{x) from plT]) . and (HSU) 



- ( E ^Ak t)A'X, + E 'I'Ak E ^Ak Bl ^Ul +T.§i^Ak t) <Ak] 
\ k k Bl k J 

nt + St 

+ Y.^Aka^Ak{^^^) / dtrT,{tr) 

ka •'^ 

= E ^Ak M\ <5t + E Bl ^Ul ^t + Y. ^Aka f 

k K Bl a 



ft+St 

dtiTaiti) 



dt' 

k 



^Ak i^^ t) SaUk +J2wAk i^^ t) (205) 



on substituting for Sa^^^^j, from Ea. (l57|) 

§-^C^Uk -Alk+J2 ^Ak Bl ^%,l + E ^Aka ^a{t+) (206) 
/ a 

and substituting from (|60p for the time derivative of the mode function. Hence 
we see that the Ito stochastic equation for ■ip'^{x,t) is 

§-^rA{x,t) -E ^Aki^.t) + E«Mfc ^l^Aki^^t) (207) 

k k 

This shows that the stochastic fields are given by 

rAi^,t)^Y'^Ukit)^Aki^^t) (208) 
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where the stochastic feature resuhs only from the phase variables being replaced 
by stochastic variables. These equations show that the stochastic field equation 
for tjj^i^ix, t) is equivalent to the stochastic phase variable equation for the a* 
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8 Appendix - Derivation of Langevin Equations 



Details in carrying out the stochastic averaging for the derivation of the Langevin 
equations in Section [2] are set out here. 
For the first order derivative terms 




E 

IJ.Ak 



-E 



da 



fj.Ak 



-F{a,a+) 



T^^Akaic 



it)) 



t+5t 



dtiTaiti) 



E 

fiAk 



da^^Ak 



F{a,a+) 



^Aki^lDn.it))^t 



(209) 



where the stochastic average rules for sums and products have been used, the 
non-correlation between the averages of functions of a^^Ak (0 time t and the 
r at later times between t to i + i5t is applied and the term involving Fq (^i ) is 
equal to zero from ([75]) . 

For the second order derivative terms, we have on expanding the product 
^'^li,Akit)^^tBiit) and using the stochastic average of a sum being the same as 
the sum of stochastic averages 



^EE 

nAk vBl 



d 



d 



da^^Ak davBi 



Fia,a+) 



2 '^t^Ak ^vBl 



d 



da^Ak da^Bl 



F{a, 



A'XkHDmmtAUc'iommt 



Bl\'^^Dm\ 



+ 



2 S/jAfc ^vBl 



Fia,a+) 



■A^ki^lDmim Eb ^BlMomit)) 1^'' dt^^t^) 



2 '^fiAk ^uBl 



-Fia,a+) 



EaKkai^lDJt))I^^''dt,r,{t,)A-i,^{al^„MSt 



_^ 2 I^Mfc T>,,Bl da^Ak dZsi 

X [Ea B'XkaKo^nm J^'' dt.T^it,) ("Id™ (^)) 1^'' dt,n{t,) 

(210) 

Using the result that the stochastic averages for the functions of the a|p„(t) 



73 



and the Ta{t+) are uncorrelated we find that 



fiAk vBl 



d d 



F{a,a+) 



KAum<Bi{t) 



liAk vBl 



d d 

do-ixAk da„Bi 



F{a,a+) 



■^Aki^lDmit))-^BlKDmit)) 



+^EE 

fiAk uBl 



d d 

da^Ak davBi 



F{a,a+) 



KkHDmimT^^BlMDmii)) 



pt+St 

X dt2rb{t2) 6t 



4ee 

nAk vBl 



d d 

da^Ak da„Bi 



F{a,a+) 



T^^AkaHDmit)) ^BlKDmit)) 



/t+St 
dtiTaih) 6t 



]ee 

/nAfe vBl 



d d 



F{a,a+) 



(■t+St ft+St 
X / dhTa{tl) / dt2Tb{t2) 



. a b 

(211) 



Now the terms involving a single F have a zero stochastic average, whilst the 
terms with two F give a stochastic average proportional to St 



/t+5t nt+St rt+St rt+St 

dtiTaih) dt2Tb{t2) = dh dt2 Ta{ti)Tb{t: 



t Jt 

/t+St i-t+St 
dti / dt2 SabS{ti - t2) 

= SahSt (212) 



so that correct to order 6t the second order derivative term is 
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juAfe vBl 



d d 
da^Ak davBi 
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= Jee 
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^EE 
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(213) 
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The remaining terms give stochastic averages correct to order St^ or higher 
that we have correct to first order in St 



F{n"(f + St). n''+(f + Sf)) - F(o^(f). a'^+it)) 



d 



da, 



IxAk 



-F{a,a+) 



nAk vBl 



d d 
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(2 
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9 Appendix - Derivation of Functional Fokker- 
Planck Equation 

This Appendix contains details in the derivation of the functional Fokker-Planck 
equation, the final form of which is set out in Section [31 

9.1 Form of Characteristic Functional 

In deriving (jlOSp we note that integrals involving cross terms between conden- 
sate and non-condensate fields are zero via mode orthogonality and hence 

dx'i'{x)E+{x) = / dxi'c{x,t)E'^{x,t)+ / dx'^Nc{x,t)E^(j{x,t) 



keC,NC 



dxE{x)¥{x) = J dxENc{x,t)¥j^c{x,t) + I dxENc{x,t)¥j^c{x,t) 



keCNC 



\ j dxEcix)E+{x)) = lY.^'^^k (215) 

•' keC 



In deriving (|113p we note that 



f dx{^Pcix,t)E+ix,t)+Ecix,t)^+ix,t)) = + ^k'xt ) 

•' ktC 

/ dxipf^c{x-.'t)'^Nci^,t) = E ^^^k I dxENcix,t)'ip+^,{x,t)^ ^ ^^.a^ 



keNC " kuNC 



9.2 Time Derivative of Characteristic Functional 

Details for the derivation of ()12ip are as follows. 

The time derivative in Eq. (|120p can be evaluated using Eqs. ()103p . (fT4|) and 
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m 

•' VktCUNC keCUNC > 

meC 

\ J dx{Y^ S,^{t)(j}^{x,t) 

7neC 



X 

rneC 



KkeCUNC keC leNC 

-^(ee^-w^^wcw|4|ee 

{keCUNC ) [.keCleNC 



ciimkmtit) (216) 

NC ) 



where the leC terms have been eliminated using Cki{t) + Cl^{t) — and orthog- 
onahty of the modes for leNC and meC has been used to ehminate the second 
and fourth contributions. 

Since for leNC and keC we have from (|103l) 



^i(t) = J dx(j)*{x,t):iNc{x,t) £,f{t)^jdx(t>i{x,t)E'^(.{x,t) 
Ut) - j dy<j,l{y,t)Ec{y.t) ^+ (t) ^ [ dy q^^{y,t)E+{y,t) (217) 



we see that 

_ld_ 
~2dt 



dx {x, t) ^ J {x,t) 



\y J '^^'^V |E E <t^*i{x.t)Ckimk{y.t)^ (zSjvc(x,t))(*S+(y,t))| 
\\j fdxdy\j2J2 (bn^,t)Ctkit)My,t)] {iEcix,t)){zE+c{y,t)) 



+ 2 



. UC keNC 



(218) 

The exponential factor in Eq. p20p can be recombined via (|107p with the 
to produce the original characteristic functional xl^c^ ^nc^ ^nc]- Writing 
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the characteristic functional in the form in Eq. (|113p and noting that multiph- 
cation by {i'ENc{x,t)){i'E.'^{y,t)) etc. can be replaced by functional differentia- 
tions such as I — \s7ri~~t) second term in Eq. (|119p is the 
sum of two contributions 



and 



/ L>2^ exp(i / ^V^^(a;,t)S"^(a;,t)) (219) 
\\j jdxdyr£Yl ^*iMCtkmk{y.t) \ {^c{x.t)){i^%c{y,t))\ 



. leC ktNC 



y" i:>2V^ exp(i J dxJ2i;''{x,t)E~^'{x,t)) (220) 



/ / dxdy \ i ^ (a;,t)C+ ^^•fcW'/'Arcfc(2/>07-+7— 



pfV-, V* 



where functional integration by parts has been applied twice. 

Equating both sides of ^x[Sc Sj, S^rc-, S^cl ^'^^ dS), (|127l) and 

P28p gives the functional Fokker-Planck equation for the hybrid distribution 
function in the form 
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10 Appendix - Derivation of Ito Stochastic Field 
Equations 

In this Appendix details for the derivation of the Ito stochastic field equations 
in Section 13] are set out. 



10.1 Derivative of Phase Space Functional Integral 

■^F[^] we first consider changes (ji^^ 



To evaluate the first term involving -S-Fl^] we first consider changes S^j^{x,t) 



in "0^4(3^1 1) during an interval St that results in a change to the functional 

From Eq. (fT82l) we have 



dt 



Hence the first term in (|I32p is given by 



Bl'^Bl 

(223) 

Using the functional Fokker-Planck equation (|129p we find after using inte- 
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gration by parts that the second term in (|132p is 

(224) 

where functional integration by parts has been used. 

10.2 Derivative of Stochastic Field Averages 

For the first order derivative terms we have on substituting from (jl39p and 
expanding 




dhTaih) 

(225) 



where the stochastic average rules for sums and products have been used, the 
non-correlation between the averages of functions of '0^(2;, i) at time t and the 

r at later times between t to t + St is applied and the term involving Ta{ti) is 
equal to zero from fTS]). Note that this term is proportional to 6t. 

For the second order derivative terms we have on substituting from ()139p 
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and expanding we find that 



fj.A,uB 



dxdy 
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dxdy 
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(226) 

where the stochastic average rules for sums and products have been used, the 
non-correlation between the averages of functions of ip'Jj^{x,t) at time t and the 
r at later times between t to t + 5t is applied. The terms involving a single T 
have a zero stochastic average, whilst from Eq. (|78p the terms with two F give a 
stochastic average proportional to St 
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(227) 
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so that correct to order 6t the second order derivative term is 



/ / dxdy ^ SrA{x,t)6rj{y,t) 
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Iff dxdy 
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\briMHBh).i)) 


\ f f dxdy 




x,y 


'w{x)N^{y)rA:B 


[srAix,t)srB'{y,t)j 



(228) 
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